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Abstract 



A theoretical and experimental investigation of the creep 
behavior of initially flat circular plates is described in 
which particular attention is paid to the effect on this 
behavior of geometry changes induced during deformation. A 
method of modeling the constitutive behavior of a material 
undergoing creep deformation in a multiaxial state of stress 
is proposed. An energy technique for solving creep problems 
using this constitutive model is described and solutions are 
obtained for the cases of clamped and simply supported circular 
plates under a uniform, time independent, lateral load. The 
technique retains the effects of both bending moments and in- 
plane forces in the deformation of the plate. A series of 
experiments ;vith lead plates is described and compared with 
the theoretical predictions obtained in this study and with 
those of previous investigators. Both the theoretical and 
experimental results show that the effect of finite deformation 
is of primary importance in the creep behavior of plates and 
must be taken into consideration even when the maximum 
deflection is well below the plate thickness. 

Thesis Supervisor: Norman Jones 

Assistant Professor of Naval Architecture 
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INTRODUCTION 



Creep is a terra coraraonly applied to a material deformation 
process in which the strains are a function of time as well as 
stress. The phenonraenon occurs in metals when they are stressed 
at high temperatures , and can result in appreciable plastic 
deformation even for stresses well below the yield stress. As 
machinery operating temperatures have risen steadily over the 
past years in search of higher efficiencies the phenomenon of 
creep has been of increasing interest to engineers. When 
designing components which will operate in an environment where 
creep may occur it is obviously necessary for the engineer to 
be able to predict creep deformations in order to produce an 
efficient and workable design. The mechanisms of creep are 
complex and as yet imperfectly understood, but enough is known 
to permit a reasonable analysis of some cases of practical 
importance. It is the purpose of this investigation to examine 
one such case, that of a thin circular plate exposed to a 
uniform lateral pressure, and to devise a means of predicting 
its creep behavior including, in particular, the influence 
upon this behavior of the geometry changes resulting from the 
deformation. It is believed that the general procedure 
developed in the solution of this particular problem can be 
applied to other structural shapes. 

Trouton and Rankine observed and described the phenomenon 
of creep in lead wire in 1904 (1). In 1910, Andrade (2) made 
the first systematic investigation of creep using wires of 
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approximately pure lead and of lead-tin alloy. In 1929, 

Norton (3) published the results of his investigation of creep 
in steel at elevated temperatures and suggested that creep 
rate, for uniaxial stress, could be represented as a function 
of stress raised to a power n. This relation, which has 
become known as Norton's law, has been the basis for much of 
the phenomenological creep analysis in recent years. Other 
functional relations for uniaxial test data have been suggested, 
such as a hyperbolic sine function proposed by Prandtl and 
endorsed by Nadai (4) and an exponential relation suggested by 
Soderberg (5). The greater mathematical tractability of 
Norton's law has made its use far more common. A generalization 
of Norton's law to the multidimensional stress state was pro- 
posed by Odqvist in 1933 (6) . Odqvist used the concept of 
octahedral stress to relate the strain rate components to the 
deviatoric stresses. At about the same time, Bailey (7) 
proposed a somewhat more complex set of relations for general- 
izing Norton's law to the multiaxial case in a form appropriate 
to the principal directions of stress and strain. Soderberg (5) 
showed how his exponential creep rule could be generalized to 
multidimensional cases through the stress and strain invariants. 
Another approach, based on the work of Drucker (8) and Hopkins 
and Prager (9) in the rigid plastic behavior of circular plates, 
uses the Tresca yield surface and the condition of normality of 
strain rate vectors to that surface to extend a uniaxial creep 
law to a multiaxial problem. This method was used by Wahl (10) 
for the analysis of rotating discs and by Venkatraman and Hodge 
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(11) for the bending of circular plates subject to creep 
deformation. At the present time there seems to be no firm 
basis, experimental or otherwise, for endorsing or rejecting 
any of the above theories as a generally applicable constitutive 
relation in creep. In this paper another model for the general- 
ization of Norton's law to multidimensional cases is proposed 

in which principal strain rates are assumed to be a function of 
the nth powers of the principal stresses. 

The analysis of the creep behavior of plates has thus far 
been confined to a bending or infinitesimal deflection analysis, 
in which in-plane forces are ignored, or a membrane analysis, 
in which moments are ignored. Odqvist (12) analyzed the bending 
of thin circular plates with clamped and simply supported edges 
using his multiaxial relations based on stress invariants. 
Venkatraman and Hodge (11) solved the same bending problems for 
circular plates using the Tresca criterion technique mentioned 
above. Venkatraman and Patel (13) (14) used this same method 
for the bending of annular plates. The Tresca technique was 
extended to the case of orthotropic bending of circular plates 
by Sandaranarayanan (15) . A membrane analysis for circular 
plates undergoing creep deformation has been done by Odqvist 

(12) (15), again using his stress invariant generalization of 
Norton's law. Onat and Yuksel (17) have suggested a very 
simple and ingenious technique for the membrane analysis of 
a thin circular plate in which it is assumed that the plate 
deforms through creep behavior into a spherical shape. The 
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stresses may then be assumed to be everywhere equal and the 
uniaxial stress-strain rate rule is applied directly. There 
has apparently been, as yet, no mathematical analysis of the 
creep behavior of plates which retains both bending moments 
and in-plate forces and which considers the influence of the 
geometry changes induced by the deformation. None of the plate 
analyses cited above were accompanied by any experimental 
verification nor have any creep tests on plates been noted by 
the author in the literature. Jones (18) (19) has shown the 
very important influence of finite deformations in the static 
and dynamic behavior of rigid-plastic plates and it would 
seem likely that creep behavior would exhibit a similar 
dependence . 

In this investigation a series of initially flat lead 
plates of three different thicknesses is subjected to various 
lateral hydrostatic pressures in a test apparatus which permits 
measurement of plate deflection at the center and at regular 
intervals along the radius. By observing deflections at 
various time intervals one is able to determine the deflection 
rate of the plate and the plate's deformation profile. Tensile 
creep tests are also conducted on specimens made from the same 
load slieets as the circular plates in order to determine basic 
material creep parameters. Experimental procedures, apparatus, 
and results are described in Chapter I. A method for mathe- 
maticallv' analyzing the creep behavior of circular plates 
including tlie effect of finite deformations is described in 
Chapter II. The method is based upon computing the rate of 
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internal energy dissipation in the plate and equating it to the ex- 
ternal work rate. It is essentially an adaptation of a technique 
suggested by Timoshenko (20) for computing the large deflection 
behavior of circular elastic plates. The problem under study 
has several inherent nonlinear aspects : nonlinear material 

constitutive relations, nonlinear geometry with finite 
deflection, and nonlinear change in creep parameters with time. 
Despite these multiple nonlinearities, by using the technique 
herein described, the results for a circular plate may be 
presented in the simple form of an algebraic equation and a 
set of graphs. All necessary computer work is incorporated in 
the graphs from which the deformation rates may be computed 
directly, with the aid of a slide rule, for plates made from 
any material whose creep characteristics are adequately 
described by Norton's law. In Chapter III the theory is 
compared with the experimental results and with the bending 
and membrane analyses of other investigators. 

Both the theoretical and experimental results demonstrate 
the strong influence of finite deflections on the creep behavior 
of circular plates. This influence is significant even for a 
deflection at the center less than the plate thickness. The 
proposed method of analysis yields results which are in quite 
reasonable agreement v/ith the experimental data. 
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CHAPTER I 



A. PLATE DEFORMATION EXPERIMENTS 
1. DESCRIPTION OF APPARATUS 

The creep behavior of a series of circular plates of 
very nearly pure lead was examined in the apparatus shown in 
Figures 1 and 2. The plates were cut from commercially 
rolled sheets of three different thicknesses commonly known 
as 4 lb., 6 lb., and 8 lb. plate. These thicknesses were: 

Plate Size Thickness 

4 lb .068" 

6 lb .105" 

8 lb .136" 

Local variations within the plates from the above thicknesses 
were no greater than ±.001". All circular plate specimens and 
tensile specimens of any one thickness were cut from the same 
lead sheet. A chemical analysis of each sheet is given in 
Appendix G. Lead was chosen for these experiments because of 
the fact that it creeps at room temperature. This fortunate 
property means that costly furnaces, temperature control 
devices, and elaborate instrumentation which would be necessary 
for elevated temperature testing are not required. Although 
lead is not a commonly used structural material we may expect 
that its behavior at room temperature will correspond to the 
behavior of more common structural materials, such as steel, 
at higher temperatures. Moreover, the theoretical analysis 
which is to be presented is valid for any material whose creep 
behavior is governed by the constitutive equation postualted. 
Thus, if agreement between theoretical predictions 
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and experiments is achieved for lead plates, then the analysis 
should remain valid for plates which are made from other 
materials whose behavior is adequately represented by the 
constitutive equation. 

The test apparatus consisted of a short section of steel 
pipe with a blank base welded to one end, and a flange and 
ring on the other end which was used to clamp the test 
specimen. The flange and ring were machined to have an 
inner diameter of exactly five inches and square edges at the 
inner lip. Centering pegs were placed in the flange so that 
the ring could always be positioned precisely above it by 
fitting the pegs into holes drilled into the ring. This was 
necessary as the inner edge of the ring must be exactly above 
the edge of the pressure vessel to insure proper end conditions. 
The flange and ring had serrations machined on the gripping 
surfaces. These serrations were relieved in both surfaces so 
as not to damage or distort the edge of the plate sample. The 
sample plate was clamped between the flange and the ring by 
eight liigh strength steel bolts tightened to a uniform tension 
with a torque wrench. Six dial indicators v;ere positioned 
above the plate with indicator holders attached to the ring. 
These dial indicators were placed so as to record the defor- 
mations at various places on the plate. Positioning of the 
indicators was facilitated by the use of extension rods of 
various lengths for the indicator contact points. All dial 
indicators were graduated to .0001" and were checked for 
accuracy with gage blocks before commencing the experiments. 
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PLATE TEST APPARATUS 

FIG. LI 
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A* Plate Test Apparatus 



B, Plate Test Sample 
(after deformation) 




C. Arransre^i ent ..f Ditl Indicat. rs 



'IG. 1.2 




A 



The apparatus was pressurized by a water column. The water 
column consisted of a length of plastic tubing connected to a 
3 diameter plastic cylinder at one end and, through a valve, 
to the pressure vessel at the other end. The use of a water 
column insured that the hydrostatic pressure in the vessel 
could be both accurately measured and accurately maintained 
over the long period of time required for a creep test. The 
cylinder on the upper end of the tube was for the purpose of 
expanding the volume of the upper part of the column so that 
small changes in the volume of the pressure vessel from plate 
deformation could not make an appreciable change in the column 
height. The pressure vessel was also equipped with a vent 
valve and a pressure gage. The pressure gage was used as a 
check to insure the pressure in the vessel was at the level it 
should be and not blocked in some way. Test pressure data 

was obtained by measuring the water column height not by reading 
the gage. 

2 • EXPERIMENTAL PROCEDURE 

The circular plate specimens, having been cut from a 
single sheet of commercially rolled lead, were flattened with 
a hydraulic press and then annealed at 200°C for two hours. 

After cooling to room temperature a sample plate would be 
mounted in the test apparatus as shown in Figure 1.3. The 
chamber of the test apparatus and the pressurizing column 
were filled with water before putting the sample in place. 

The sample plate was marked by pencil with concentric circles 
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from the center at 1/2 inch intervals, and then placed 
carefully on the upper flange of the test apparatus. The 
securing ring was then put into position above the specimen 
and secured with eight bolts tightened to a uniform tension 
with a torque wrench. The dial indicators, having been 
previously attached to the securing ring, were then positioned 
to their desired locations at the center and at intervals 
along the radius. The assumption was made that deformation 
was axisymmetric so that readings would be a function of 
radius only. Observation of the tests revealed no evidence 
that this was invalid. With the water column adjusted to a 
height giving the desired hydrostatic pressure the test is 
ready to begin. Before releasing the pressure, initial 
readings are taken on all dial indicators. These are the flat 
plate, or t=0, readings to which all subsequent readings are 
referred. The test is started by opening the valve between 
the pressure column and the pressure vessel, thus admitting 
a hydrostatic pressure to the back of the plate. Dial 
indicator readings and room temperature are subsequently 
recorded at various time intervals. Plotting these readings 
at any given time gives the deformation profile at that time, 
and plotting the readings at any given radius over all times 
gives the deformation history. The slope of the deformation 
history at any particular time is the deformation rate at that 
time. Of particular interest is the largest deformation which 
occurs at the center. Two tests were made at each pressure 
using a different plate of the same thickness in order to check 
the reproducibility of the results. 
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3. 



RESULTS 



The results of the plate deformation experiments are 
summarized in Figures 1.4 through 1.16. Raw data is listed 
in Appendix F. 

Figures 1.4 through 1.6 are plots of the deflection of 
the center of the plate as a function of time. Two ordinate 
scales are shown, one for the deflection in inches and 
another for deflection as a ratio of the corresponding plate 
thickness . 

Figures 1.7 through 1.9 also show center deflection as a 
function of time, but here all deflections are plotted relative 
to that at t=l hour in order to show more clearly the relative 
deflection rates for different pressures. The time t=l hour 
was chosen arbitrarily as a reference point. Some other point 
could have served equally well. 

Figures 1.10 through 1.16 show representative deflection 

w 



profiles for increasing values of 
verse deflection as a function of 
time t=l hour. Similar plots for 
show no significant differences, 
plots that the deflection profile 
curves of the forms 

w = w (1-^)^ 
o a 



The plots are of trans- 
radius . Those shown are at 
other times such as t=5 hours 
It can be seen from these 
is approximately bounded by 

1.1 



and 



w = w (1 — r) 
o a^ 



1.2 



where w^ is deflection at the plate center. For tests at the 



w 



lower 



ratios, the measured profiles correspond quite closely 
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Relative plot of 8 lb. plate results 
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w 

to the former expression, and at higher ratios the observed 
profile seems to migrate toward the latter parabolic form. 

This seems rather reasonable as the former expression corresponds 
to the profile given by Timoshenko (21) for the elastic 
deformation of a clamped, hydrostatically loaded, circular 
plate and the latter is the corresponding membrane profile. 

The fact that the plate behaves in this way was of particular 
significance in suggesting the method of analysis described in 
the next chapter. 



B. TENSILE CREEP TESTS 
1. DESCRIPTION OF APPARATUS 

Tensile creep tests were necessary for obtaining the 
appropriate creep constants, that is, K and n in Norton's law 
v/hich may be written as; 

e = Ka^ 1.3 

Tensile test specimens were prepared from the same lead sheets 
from which the plate specimens were made. The tensile 
specimens were also flattened with the same hydraulic press 
and annealed for two hours at 200°C in a manner similar to 
the circular plates. The type of specimen used is illustrated 
in Figure 1.17D. 

The tensile testing apparatus is shown in Figures 1.17 
A, B, C, and E. It consists essentially of a frame in which 
the weighted tensile specimen is suspended and a cathetometer , 
an optical measuring device, for measuring the distance 
between two marks inscribed on the specimen. There is also 
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a modified hydraulic jack used to facilitate the rather 
important task of smoothly loading the specimen at the 
commencement of a test. The specimens were loaded with 
calibrated iron weights mounted on a pan suspended from the 
lower specimen's grip. 

The frame was constructed of aluminum I beams bolted 
together. The specimen was suspended from an eye mounted 
in a transverse bar centered in the top of the frame. A 
flexible braided wire with an eye soldered to each end was 
used to suspend the specimen. Each of these eyes had its 
inner edge machined to a knife edge to insure that a specimen 
would hang freely. A slotted bar, as shown in Figure 1.17C, 
was used to prevent any rotation of the specimen during the 
test. This was necessary to insure the marks were always 
visible by the cathetometer . The bar was placed loosely on 
its supports and provided only a very gentle, yet sufficient, 
restraint against any tendency for the specimen to rotate. 

The specimen was held between two grips as shown in 
Figure 1.17C. The upper grip was made of aluminum and the 
lower of plastic. The use of plastic was to minimize the 
weight of the lower grip. This assisted in preventing any 
bending of the specimens while mounting them in the frame. 

The initial loading of the specimens was accomplished 
with tlie aid of a modified hydraulic jack. Prior to the 
commeiicement of a tost the weights and weight pan were 
supported by a small hydraulic jack as shown in Figure 1.17E. 
The jack is adjusted so tliat the hook on the weight pan is 
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barely out of contact with the ring on the lower grip. When 
the jack's bypass valve is released, two springs draw the 
piston quickly and smoothly down from the weight pan, thus 
loading the specimen. As seen in Figure 1.17B the radial 
slots in the weights ride on a verticle rod. This rod 
prevents the weights from rotating and twisting the specimen 
but does not inhibit vertical movement of the weights. 

Measurement of tensile strain was accomplished with the 
aid of the optical device shown in Figure 1.17 A, B, and C, 
known as a cathetometer . It consists of two telescopes 
mounted so that they can be moved vertically as a unit or 
moved apart in the vertical direction. The unit had a vernier 
for measuring a change in the vertical position of the two 
telescopes as a unit. Each telescope was equipped with a 
movable hairline whose deflection from its centered position 
was measured by a vernier graduated to .0001". By following 
the separation of two scribed marks initially a known distance 
apart on the specimen, one can use this instrument to measure 
creep elongation and have no physical contact with the 
specimen which could possibly affect the normal deformation 
process . 

2 . EXPERIMENTAL PROCEDURE 

The tensile specimens, having been prepared from the lead 
sheet as described in the previous section, were scribed with 
two thin max'ks one inch apart on the narrow section of the 
sample. The scribe marks were made by clamping a one inch 
gage block over the specimen and drawing a razor blade down 
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either side of the gage block. Great care was necessary in 
attaching the grips to the specimen in order to avoid bending 
the soft lead. A specially prepared block, machined to the 
proper dimensions, was found useful in this regard for 
supporting the specimen and grips while they were being 
attached. The slotted bar, used to prevent rotation of the 
sample, was placed around the specimen when attaching the 
grips. This assemblage is then suspended in the frame by 
placing the hook on the upper grip into the eye at the lower 
end of the braided wire. The slotted bar rests on ledges 
placed on either side of the frame. The weight pan is attached 
to the lower grip and the jack is raised to support the pan. 

The required weights are then put in place on the pan. The 
telescopes are adjusted so that, with the crosshairs centered, 
the crosshair in the upper telescope is aligned with the upper 
scribe mark and that of the lower telescope with the lower 
scribe mark. The test is begun by releasing the jack and 
allowing the weight to be applied to the lower grip. As the 
specimen strains the scribe marlcs move a\;ay from the cross- 
hairs. The amount of elongation can be measured by keeping 
the centered crosshair in the upper telescope aligned with 
the upper scribe mark and moving the lower telescope crosshair 
to realign it witli the lower scribe mark. The amount the 
lower crosshair is moved can be read from the vernier. The 
scribe marks, although quite thin to the naked eye, are much 
wider than the hairline when viewed through the telescope. 

The accuracy of measurement must therefore be considered as 
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about 1.0002" rather than within .0001" which is the 
graduation of the vernier. 

By measuring the movement of the scribe marks over a 
period of time the deformation history of the sample under 
the given load can be determined. Performing a series of 
these tests at different loadings yields a family of curves 
showing the deformation history at different stress levels. 

By measuring the slope of each curve at some given time one 
can obtain a strain rate, for that time, at each stress level 
tested. Plotting strain rate versus stress on a log-log 
plot will, in the case of a material whose creep properties 
conform to Norton's law, yeild a straight line from which the 
creep constants K and n can be determined. 

3. RESULTS 

Strain versus time for the 8 lb., 6 lb. and 4 lb. plates 
is plotted in Figures 1.18 and 1.19, Figure 1.20, and in 
Figures 1.21 and 1.22 respectively. The slopes of these 
plots are measured at times t=l hour and t=5 hours. These 
two times wore selected as t=l hour is a point apparently 
well within the primary creep region where deformation rate 
is changing rapidly; and, by constrast, time t=5 hours has 
an appreciably smaller and nearly constant creep rate more 
characteristic of secondary creep. Creep rate versus stress 
is plotted on a Jog-log scale in Figures 1.23 through 1.25. 
These data indicate appropriate creep constants as follows: 
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8 lb. plate n=3.7 



K = 5.56 X 10 @ t=l hr. 

K = 2.30 X lO"^"^ 0 t=5 hrs 



6 lb. plate n=5.0 



— 1 9 

K = 4.10 X 10 0 t=l hr. 

_ 1 O 

K = 1.70 X 10 0 t=5 hrs 



4 lb. Dlate n=4.0 



K = 8.30 X 10 0 t=l hr. 

K = 3.50 X lO"^^ 0 t=5 hrs 



It is interesting to note the marked differences in the 
creep characteristics of the samples cut from the three 
different sheets. All three sheets were made from lead of the 
same commercial grade and one might have expected them to have 
exhibited nearly identical creep characteristics. Such was 
not the case, however, and the differences are too great to 
be accounted for by experimental error or data scatter. 

Gifkins (22) made a careful study of the effect of alloying 
on the creep of lead using lead- thallium alloy. He found 
that, in very dilute alloys, a change of composition of as 
little as .01% thallium would result in a significant change 
in grain size and creep performance. Creep characteristics, 
therefore, certainly in the case of lead alloys and probably 
also for other materials, seem to be quite sensitive to small 
changes in the composition of the material. This means that 
calculations based on assumed creep constants may be greatly 
in error unless the constants have been derived from tests on 
exactly the same material as that being considered. 
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TIME (HRS) 
FIG. 1.19 
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TIME-(HRS) 
FIG. 1.21 
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FIG. 1.22 
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8 LB. PLATE 




Log-'Log plot of Tonsi-le Test Results, 
8 lb. rlate 
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6 LB PLATE 




FIG. 1.24 

L'^g-Log plot of Tensile Test Results, 
6 lb. plate 
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4 LB. PLATE 




Lf^'r-Loc "lot of Tensile Test Results, 
4 lb. plate 
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CHAPTER II 



Theoretical Analysis 



A. ASSUMPTIONS 

Consider the case of a homogeneous, incompressible, 

isotropic, circular plate subjected to a time-independent, 

lateral, hydrostatic pressure. The plate is assumed to be 

sufficiently thin for the problem to be treated as plane 

stress with a =0. We assume the creep behavior of the 
z 

material may be modeled in the manner discussed in Appendix A, 
which for this case gives: 

= K[|a^|^ sgn - -^lagj^sgn Og] 2.1 

= K[|ag|^ sgn Cg - sgn o^] 2.2 



where sgn is the signum function defined as 



+1 x>0 
^ = [-1 x<0 

In this form 2.1 and 2.2 can be used for negative as well as 
positive stresses. 



B. BASIC EXPRESSIONS FOR MOMENTS AND IN-PLANE FORCES 

Solving 2.1 and 2.2 for and Og in terms of the strain 
rates we find: 



(— ) 




1/n 


sgn 


(Jr+i^0> 


2.3 


^ 1/n 




|l/n 


sgn 




2.4 



The extensional strains along the middle surface of the 
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1 



6 and the curvature of the 



plate will be denoted as e^ and e 

neutral surface as and Kq. Adopting the Kirchoff assumption 
that sections which are plane before deformation remain plane 
during deformation, the total strains, and e^, at any point 
along the thickness of the plate may be written as: 



e 


= e 


+ 


ZK 


2.5 


r 


r 




r 




^0 


II 

CD 

CD 


+ 


ZKg 


2.6 



where z is the distance from the middle surface. Notation and 
sign conventions are as shown in Figure 2.1. Taking derivatives 
of 2.5 and 2.6 with respect to time gives the strain rates: 



e = e + ZK 
r r r 



=6 ' ®e + ^*6 



Substituting into 2. 
„ /4 . 1/n I . 

Or = (3K> hr 



_ _ ,4 s 1/n I 

°e I’ 



3 and 2.4: 

+ + Z (K^ + i'Kg) I 

sgnte^+ieg +z ( k^+^Kq) ] 

+ -ze^ + z (Kq+^Kj.) 1 

sgn [Og+le^+z (Kq+^Kj.) ] 



2.7a 



2.7b 



2.8 



2.9 



From this point conventional notation for the signum function 
will become rather onerous and we shall define the symbol 

Xj, x^) |P 

as meaning the absolute value of the function f is raised to 

the power p, the result is then assigned the sign which 

function f (x, , x„ , x ) assumes when evaluated with the 

12 n 
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appropriate arguments , x^, x^) for the case in 

question . 

Considering a homogeneous thin plate we may compute 
in-plane equivalent forces per unit length, and N^, as: 



N 



N, 



f- 

2 a^dz 



2 

. 



h 

2 a 
h 
2 



dz 



2.10 

2.11 



Equivalent moments per unit length are: 



M 



M, 



h 

’2 

2 

rh 

2 

h 

~2 



a zdz 
r 

a a zdz 

V 



2.12 



2.13 



By substituting (2.8) and (2.9) into the above equations, 

(2.10) through (2.13), and integrating over z, we may obtain 

expressions for the plate force and moment resultants, , 

N„, M , and M-. Some care, however, must be exercised in 
Ur u 

performing this integration. The signs of various terms in 
the result depend upon the relative magnitude and directional 
sense of the extensional strain rate and the curvature rates 
in the plate. It can be shown, however, that upon integration 
all cases may be represented by expressions of the form: 
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1 /^ _ JL _ 
'3K^ n +1 



K + 7IC 

r 0 









2.14 



A.l/n n 
^3K^ n +1 



h/.^ ^i.-. ^ |l+l/n 



K ^+ 7 K 

0 ^ r 



., 0®e 2<'^e^K 






2.15 



,4 . 1/n n 1 rh I . . 1 . .h,. ,i. . 1 1+1/n 

(^) — ; [9 |e +ie +^(K +iK.) I 



■ 3 K ' 



n +1 . . 1 . ‘■2 ' ''r ' '^^ 0 ' 2 '■'r ' ^■'0 

■=r 



n 



( 2 n+l) (Kj.+iKQ) 0 2 r 0 s 



^ ^ I 2 + 1 /n 






n 



( 2 n+l) (k^+^Kq) 






2.16 



/4_x 1 /n _n 

'3K^ n +1 . 






•||®9 



n 



( 2 n+l) (iCp+^Kj,) 



l-e 






n 



( 2 n+l) ('^ 0 +i'^j-) 



^ + 1 ^ _ li(|^ +4.,^ v| 2 +l/n, 

0 ^ r 2 '^0 ^ r^ I s ^ 



2.17 



c. DEFORf'lATION PROFILE AND STRAIN RATES 



For a circular plate undergoing finite deformation the 

curvatures and extensional strains can be written: 

w ' 



K = -W 

r 



0 r 



2.18 



e = u* + 



(w' ) 



> \ 2 



u 

®0 ■ r 



2.19 



where w and u are the transverse and radial deflections 
respectively with directional sense as shown in Figure 2.1. 
Taking derivatives with respect to time we obtain the rates: 



IC = "W 

r 



. w ' 

“ “r 



2.20 



* * I ■ * I I 

e = u + w' w 
r 



. u 

®0 “ r 



2.21 



To proceed further, the transverse and radial deformation 
rate profiles for a plate with a given set of edge conditions 
must be determined. For the finite linear elastic deflections 
of a clamped circular plate Timoshenko and Woinowsky-Kreiger (23) 
employed the elastic bending analysis profile: 



w = w (1-^) ^ 



2.22 



where w is the deflection at the center and a is the outer 
o 

radius. In the experiments described in Chapter I it was 

cl 

noted, for the r - ratios examined, that the measured deformation 
h 

profile could be represented by equation (2.22), especially at 

w 

lower values of Assuming, therefore, that the vertical 

deformation profile may be represented by (2. 22) the curvatures 



and curvature rates may be written: 



4w 






2 ( 1 -^-) 



4w 



'^0 






2.23 
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2.24 



K = 

r 



4w ,2 



4w 2 

^6 = -;?<i-|-> 



For the deflection rate in the radial direction, u, we 
adopt an expression suggested by Timoshenko and Woinowsky- 
Kreiger (24) for large deflections of clamped circular 
elastic plates which, adapted for time dependent problems, is: 

2.25 



u = r(a-r) (c, + c„r + c_r^ + c 

12 3 n-i 



This expression satisfies the end conditions u=0 at r=0 and 
r=a and provides for a wide variety of possible profiles 
through variation of the c^ values. The constants, c^, will 
be determined by requiring minimum internal energy dissipation 
in the plate. Taking the first two terms of the series we 
have : 

O O'! 

2.26 



2 ^ 2 3 

u = Cj^ar - Cj^r + 



from which: 



u' = (a-23:)cj^ + (2ar-3r )C 2 



u , . ,2. 

- = (a-r)Cj^ + (ar-r )C 2 



2.27 



2.28 



With these expressions the extensional strain rates, (2.21), 
become : 



e = (a-2r)c, + (2ar-3r )c„ + 
r 1 2 



2. _ . ,, r^2 



d 



2.29 



e^ = (a-r)Cj^ + (ar-r )C 2 



2.30 



Examining the expressions for e^ and e^ we are led by 
dimensional considerations to suspect that Cj^ and C 2 are of 
the form: 
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2.31 



c 



1 



w w 
o o 



w w 




where c| and are pure numbers. In order for each term in 

the sum to have the time dimension each must contain . If 

e^ and eg are to be zero at any point on the plate when 

w =0, then each term should contain w . To be consistent 
o o 

with the last term in (2.29) the diameter, a, is selected as 

a length factor. If one considers extensional energy 

dissipation only and solves analytically the n=l case for c^ 

and c^r it is seen that they are indeed of the above form. 

_ 

Defining x = — , we may now write the curvature rates and 

extensional strain rates as: 

4w 

K = — ^ (1-3x2) 2.32 

r 




(1-X2) 



2.33 



e 



e 



r 



9 



w w 
o o 

-"T-2 — 



W W 

o o 

—T2 — 



[(l-2x)c^ + (2 x-3x2)c 2 + 
[(l-x)cj^ + (x-x 2 )c 2 ] 



16 x 2 (l-x 2 ) 2 ] 



2.34 



2.35 



We now have expressions for the curvature rates and exten- 
sional strain rates in terms of geometrical parameters for a 
clamped circular plate. These expressions can now be sub- 
stituted into 2.14 through 2.17, the force and moment resultants. 

w 

For the higher values of the experiments showed that the 

deformation profile tends to the parabolic membrane form: 
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The problem may be solved using this profile also and, as 

will be seen, the solution for the center deflection rate, 

w , tends to converge to the same value for the two profiles 
° w 

as increases. The same procedure can be done with other 
assumed profiles as seen in Appendix D for the simply 
supported case. 



D. ENERGY EXPRESSIONS 



The total internal energy dissipation in the plate during 
creep deformation may be expressed as: 

■2tt ra 



V = 



(M^K^+M^i<:g+N^e^+NQeg)rdrd9 



2.37 



'0 0 

where transverse shear energy dissipation is neglected. The 
integration, strictly speaking, extends over the deformed 
surface area of a plate. However, if attention is confined to 
moderate transverse deflections and small strains then the 
deformed surface can be taken to equal the original surface 
area. This procedure is usually adopted for deriving a 
consistent set of equilibrium equations and strain relations 
when using the principle of virtual work. (25) Using the 
previously derived relations, equations 2.14 through 2.17 and 
2.32 tlirough 2.35, we may now write for the clamped case: 



56 



,, . ,4 X 1/n n ,2+1/n ,'^o. 1+1/11,1 i „ . r> l 

M K = (t^) —nr h ' (— Z-) B+A 

r r 3K n+1 a ' ' 



1+1/n 



n 



(2n+l) A 



7 I B+A I i |b-a 1 



n 



(2n+l) A^ 



b-a| 



2+1/n 



] ' [l-3x^] 



2.38 



M K = _n_ ,2+1/n 1+1/n 1. > 1+1/n 

^09 ^ 3 k ’ n+1 " 



(2n+l)D^ |E+D|s + d'^ 



n 



(2n+l)D' 



|Ef-D|2+l/nj . 



2.39 



M \ 1/n n ,2+1/n , o, 1+1/n 1 1+1/n 

Vr ' <3K> STl'’ 2A 

-|b-a| 1*1''"1 • [(l-2x)cj^ '^ + (2 x-3xMc^ ^ 



w 

+ 16 x^ ( 1 -x^ ) ^ ] 



2.40 



. .4 V 1/n n ,2+1/n /'^o. 1+1/n 1 r i t^. . r> I 



2D 



w 



* [(I-X)c^ ^ + (x-x 2 )c^ 



2.41 



where : 
A 
B 



(3-7x2) 

w w 

(i - ix)c' ^ + (|x - lx2)c^ ^ + 16x2(1-x2)2 % 



D = ( 3 - 5 x 2 ) 

E 



w w w 

= (j - 2x)c| + (2x - 1x2)0^ + 8x2(l-x2)2 _ 
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with these expressions 2.37 can be written as: 



• _ 2/4 V 1/n n ,2+1/n ,^Ovl+l/n 

V = 27raM3K) " HTT h 

/I (f,, + f,, + f„ + )xdx 2.42 

0 ^4 M- N N- 

y -r A -r A 



where the quantities f are all functions of x, c,' , c', 

3., -L 2- 

b 

w 

and n. These symbols, f , represent the terms containing 

b 

X in 2.38 through 2.41 which must be integrated over the radius. 

The result of this integration over x will be a function of 
w 

c|, c^f j^/ and n which will be denoted I. The equation for 
V may now be written as : 



V = 2na' 



(— ) 



n 



n+1 



2+1/n 



,'^Ov 1+1/n 



w 



I (c| ,c 



o 



2'h 



,n) 2.43 



The quantities c| and c^ as previously discussed determine 

the radial deformation rate profile. They must be chosen such 

that the rate of energy dissipation in the plate is minimized. 

Clearly, minimizing V with respect to cj^ and c^ is equivalent 

to minimizing I with respect to c^ and cl. The minimization 

w 

of I must be done for the current value of ^ and n. The 

h 

correct values of cj^ and c' for one set of values for the 
w 

parameters ^ and n may, in general, be different from those 

for some other set. One way of determining the appropriate 

w 

values of c^ and c^ might be to substitute values for and n 
into an expression for I and take derivatives with respect to 
and c^, set them equal to zero, and solve for the values 
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resulting in a minimum I. Unfortunately, the expression 



is such that the integral cannot be evaluated analytically 
so there is no analytical expression for I which can be 
differentiated with respect to c| and c^- We must, therefore, 
evaluate I numerically and search for those values of cj^ and 
c '2 which result in the miminum value of I for a given case 
(w^/h, n) . This was done, for the cases considered, with the 
computer programs whose listings are given in Appendix H. 

These programs were used with an IBM 1130 Computing System 
at Massachusetts Institute of Technology. The computer 
output yielded a series of smooth curves for regular variations 
of c| and c^. These curves showed clearly identifiable 
minima with respect to cj^ and c^ for the function I. The 
procedure is illustrated for a typical case in Appendix E. 



E. SOLUTION FOR DEFORMATION RATE 

If a structure is assumed to deform in a manner prescribed 
by a kinematically admissible virtual velocity field, (i.e., 
continuous and satisfies the velocity boundary conditions) , 
then the principle of virtual velocities states that the total 
virtual internal energy dissipation is equal to the total 
virtual external energy dissipation. 
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n) we 



w 

Having minimized the integral I for some case (j^/ 
find the deformation rate, w^, by using the principle of 
virtual velocities and considering a small variation of the 
deformation rate about its equlibrium value. This gives; 



^ 6 w^ = 2tj f^q6w rdr 



2.44 



where q is a time-independent, transverse pressure. For the 
case of the clamped plate we have: 



2ir /^q6w rdr = 2irq6w /^(l-^T)^rdr = ■ ■ 6w 

0 O 0 ^ Jo 



2.45 



and : 



dV I, . . 1/n 2/4 .l/n,l . 1+1/n, 2+1/nT i x ♦ 

■ 3 -!- 6w - w ' [ 2 Tia^{- 3 rr) ' (—t) h I] 6 v; 

dw o o 3K a'^ o 

o 



2.46 



The expression 2.46 is the derivative of V for some specific 
w 

value of and the derivative is properly taken considering 
w 

as a constant. We are seeking the value of w^ for a given 
w 



instant when 



has some fixed value. An instant later w 



h ■ ' ' ■ ■' ■■ ' o 

will, of course, have some other value and at that time w^ 

will also have some other value, but at a given instant w^/h 
may be considered fixed. Equating 2.45 and 2.46 we obtain 
the result: 






2 + 2 /n 



1,2+1/n „ 



n 



2.47 



or 



w 



r,a.2+l/n ,3Ka.l/n ,1 ..n 

= [ (u ) (-T-) q (^) ] 



'h 
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2.48 
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If we define the quantity $ as: 



w 



$ = [ 



( 



o 

Kaa>^ 



) 



3/4(3)2n+l 



1/n 



2.49 



This can be written in the non-dimensional form; 



" 3/4 (jj) 4 



2.50 



The parameter is a direct function of I, the form of this 
function depending upon the assumed velocity profile. For the 
clamped case: 



$ 




2.51 



It is convenient to define 4> as it permits a direct correlation 

of results among the different velocity profiles considered. 

When the quantitv I has been minimized with respect to c| and 

, w 

c^f is a function of n and ~. Values of <I> , for the clamped 

w 

case, are plotted versus various values of n in 

Figure 2.2. Use of these values with any of the equations 2.47 
through 2.50 will give the center deflection rate w^ as a 
function of the applied pressure and the plate's material and 
geometrical oroperties. Once w has been determined, the 

.. i. Q 

deflection rate at any other point on the plate can be found 
from the assumed velocity profile, which in this case is given 
by equation 2.22. 
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FIG. 2.2 



It was noted in the plate experiments that for the 
w 

higher values of the profile began to tend toward a 

parabolic form (or membrane profile) described by the 

equation 2.36. The entire analysis just described can be 

repeated using the profile 2.36 vice 2.22 to obtain values 

of ^ for the parabolic configuration. When this is done it 

is found, as shown in Figure 2.3, that the difference between 

w 

the two solutions is small for pp- = 1.5 and has virtually 

w 

disappeared v/hen = 3.0. It seems, therefore, that when 

w 

the plate deforms to a value of p^ large enough for there to 
be a significant deviation from the profile of equation 2.22 
it will have deformed enough so that the solution for w^ is 
relatively insensitive to the assumed profile. 

The computer program listed in Appendix H is written for 
an assumed velocity profile of the form: 

w = * [a, (1-^)2 + a,(l-|^)] ' 2.52 

where + a 2 = !• The analysis utilizing this profile is 

outlined in Appendix C. If we assign the values = 1 and 
= 0, the analysis with 2.52 is the same as for the clamped 
prof ile, 2 . 22 . If we assign values = 0 and «2 = 1 the 
analysis is for the parabolic membrane profile, 2.36. Some 
experimentation with fractional values for and was 

attempted, but these produced no more useful results. The best 
aqreement with experimental data came with and at their 

extreme values as above. The final result using this profile 
is : 
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n = 4.0 



.20 i=. 



• • a ,= 1.0 CLAMPED PROFILE 

X X a, = 0.0 MEMBRANE PROFILE 



.15 



.10 



.05 




0 1.0 2.0 3.0 4.0 

h 

FIG. 2.3 

Co'’it>arison of clarnoed and rrenbrane solutions 
for n = 4.0 
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Kaq c 



2.53 



so in this case 
3-a, 



$ = ( 



12 I 



-) 



2.54 



where is the appropriate integral for the combined profile. 



F. SIMPLY SUPPORTED CASE 

The simply supported circular plate may be treated in 
the same manner as above if the velocity profile appropriate 
to it is known. In the elastic bending analysis of this case 
Timoshenko and Woinowsky-Krieger [26] determined a deflection 
profile of the form: 



^ 1+v r^ 



W — W (1 2-) (1— p 

o a^ 5+v a 



) 



2.55 



which for v = 1/2 gives the velocity profile: 



r ^ 3 r ^ 

w = w (1-^) . 

o a'^ 11 a“^ 



2.56 



There is no experimental data available in order to verify 
this profile as in the clamped case, but at least a reasonable 
estimate of the creep behavior of the simply supported case 
can be expected from an analysis with the profile of 2.56. 

This analysis is outlined in Appendix D and yields a final 
result of: 



w 



Kaq 



-) = 3/4 (^) 



n 



‘22 I 



2.57 
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so that in this case: 

t = (jA-) 2.58 

s 

where I is the integral appropriate to the simple support 

w 

profile. Curves of <I> versus ^ for various values of n are 

h 

given in Figure 2.4. In Figure 2.5 the results of the 
clamped and simply supported cases are compared for the 
case n = 4 . 
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CHAPTER III 



Comparison of Theoretical and Experimental Results and Results 

of Previous Investigations 

A comparison of the experimental and theoretical results 

is shown in Figures 3.1, 3.2, and 3.3. Here the non-dimensional 
w 2n+l 

quantity ( r\^ ^ plotted versus the non-dimensionalized 

Kaq w 

center deflection, for the creep parameter n appropriate 

to each of the three types of plate tested. In these plots 

the appropriate analytical results for comparison with the 

experiments are, of course, those for the clamped case. Also 

shown are plots of the bending analysis results of Venkatraiman 

and Hodge [11] , and of Odqvist [12] ; and the membrane analyses 

of Odqvist [12] [16] , and of Onat and Yii ksel [17] . The 

mathematical details of reducing these various analyses to a 

form appropriate to this plot are given in Appendix B. The 

bending analysis solutions are horizontal lines since that 

type of analysis does not consider geometry changes. The 

membrane analyses consider geometry changes but not bending 

moments within the plate. Experimental points are given at 

each pressure tested for times of 1 hour and 5 hours. These 

points are found by obtaining the deflection rate at the 

indicated times from the plate test experimental data and 

substituting this value along with the appropriate creep 

constants found in the tensile tests and the pressure and 

geometrical constants into the form ( --) / (^) ^ . The 

Kaq ^ 
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2n*\ 




0 .5 1.0 1.5 



h 

FIG. 3.1 

Comparison of experiment and theory , 
8 lb. plate 
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2n +1 



VENKATRAMAN a HODGE 



ODQVIST BENDING 




Wo 

h 

FIG. 3.2 

Comparison of experiment and theory, 
6 lb. plate 
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I +u^ 




h 

FIG. 3.3 



Comparison of experiment and theory, 
4 lb. plate 
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value thus obtained is then plotted versus the appropriate 
w 

value of as determined from the plate test data. The 
times t=l hour and t=5 hours were chosen as representing 
cases of rather rapid and relatively slov/er deformation 
respectively. 

The marked influence of finite deformation on the creep 
behavior of a clamped circular plate is clearly indicated by 
both the theoretical prediction of equation 2.50 and the ex- 
perimental results, which are seen to be in substantial 
agreement. 

A similar plot, Figure 3.4, shows the relationship of 
the simply supported case analysis to the results of the 
simple support bending analysis of Venkatraman and Hodge [11] 
and of Odqvist [12] and the membrane solutions of Odqvist [12] 
[16] and of Onat and Yii ksel [17] , for the case n = 4. No 
experimental data is available for this type of end condition, 
but the analytical results also indicate a marked influence 
from finite deformation. 
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Wo 

h 

FIG. 3.4 

Comparison of simply sunrortod olate theories 
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Summary and Conclusions 



Experimental and theoretical results have been presented 
which show that the effect of finite deformation upon the creep 
behavior of circular plates is of primary importance. It is 
necessary to consider this effect even when the maximum 
deflection of the plate might be well below the plate thickness. 

A method of modeling the constitutive behavior of a 
material undergoing creep deformation under a multiaxial state 
of stress has been proposed. A technique for solving multi- 
axial creep problems by using an energy method and the 
suggested constitutive equations yields a result in the form 
of an equation and graphs which can be quite easily applied. 
Solutions were obtained for a circular plate under a time 
independent uniform lateral load with clamped and simply 
supported edge conditions. A comparison of theoretical and 
experimental results for the clamped case showed substantial 
agreement; indicating that the analytical technique described 
may be expected to provide quite reasonable estimates of plate 
deformation rate when the creep behavior of the material may 
be adequately represented by Norton's law, and the appropriate 
values of the material creep parameters are known. A comparison 
of these results with those of several previous investigators 
has also been presented. Lead was selected for the experimental 
material in these tests but the analytical results should 
remain valid for any material whose behavior is adequately 
represented by the constitutive equations. 
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Specific consideration was given only to circular plates; 
it can be expected, however, that the conclusion as to the 
importance of finite deformation will apply as well to plates 
of other shapes, initially straight beams, and to initially 
straight cylindrical shells. 



76 



APPENDIX A 



Constitutive Equations for Multiaxial Creep 



Creep data is not plentiful, and that which is available 
is almost always the result of uniaxial tests. It has been 
found that much of this experimental creep data can be repre- 
sented by an expression of the form: 

e = K (A.l) 

where K and n are material parameters which are, in general, 
functions of time and temperature. This expression is 
commonly known as Norton's law. [3] 

In generalizing Norton's law to a form suitable for 
multidimensional problems the most common approach has been 
to use the concept of equivalent stress and equivalent strain 
rate for determining the proportionality between strain rate 
and deviatoric stress. Mendelson [27] uses a relation of the 
type : 



2 



1 ] 



(A. 2) 



where creep equivalent strain rate, ^ , may be related to the 

e 

equivalent stress, a , by Norton's law. 

e 



Mere : 



Op = /3J2 (A. 3) 

* 

and is the strain rate resulting from a stress of this 
magnitude. Odqvist [6] uses the expression: 

“i] 2 c 'a ' o (A. 4) 

c c 

where is some standard value of creep rate for which 0 = 0 ^. 
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These expressions seem reasonable and satisfy the mathematical 
requirements for raultiaxial creep constitutive relations , but 
are often rather intractible for analytical work since one 
must always compute in order to know the proportionality 
between the local stress and creep rate. 

In this analysis the constitutive behavior of a material 
undergoing creep deformation in a multiaxial state of stress 
will be modeled in the following way; 

^i cFj -V (A. 5) 

where the subscripts i, j, and k represent principal directions 
of stress and strain rate. It can be shown that for conser- 
vation of volume we must have v = 1/2. In real materials such 
things as voids and inclusions may give rise to situations in 
which volume is not conserved. If this effect is important 
then in such cases some other value for v must be assigned as 
appropriate to the specific material under consideration. 

It is also desirable that the constitutive relation (A. 5) 
hold for both positive and negative a for all values of n. 

The significance of positive and negative signs here is merely 
one of direction, tensile or compressive, and this directional 
sense must be retained for all values of n. This may be 
accomplished through the use of the signum function defined as; 



r +1 X > 0 

=9n X = { _i ^ ^ g 



(A. 6) 



An expression of the form Ix|*^ sgn x is used to denote that the 
absolute value of x is raised to the n power and the result is 
assigned the sign of x. 
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For homogeneous isotropic materials, therefore, the 
proposed constitutive relation takes the form: 

sgn sgn sgn (A. 7) 

These equations are valid for a material which has the 
following properties: 

(1) Homogeneous, incompressible, isotropic material. 

(2) Coaxial stress and strain rate tensors. 

(3) Zero deformation rate under hydrostatic pressure. 

(4) Coordinates correspond to principal directions of 
stress and strain. 

(5) Creep behavior corresponds to Norton's law in 
uniaxial stress. 

It should be noted that this model predicts an effect from 
a hydrostatic pressure superimposed on a deviatoric stress. 

The validity of this may be questionable and was not really 
tested in the cases studied in this investigation. For 
applications where larger hydrostatic pressures are involved it 
may be better to use only the deviatoric components of the 
stresses . 

For the case of a thin circular plate, choosing cylindrical 
coordinates and assuming = cr =0, equations A. 7 become: 

£r = ^r *^0^ 

Eg = K[|ag|" sgn sgn a^] (A. 8) 

Eg " sgn sgn Ug]. 
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APPENDIX B 



Interpretation of Results of Previous Investigators 
Introduction 

The following theoretical treatments all apply to circular 
plates loaded with a time independent pressure uniformly 
distributed over the entire plate area. The infinitesimal 
solutions in Sections I and II retain only bending moments 
while solutions in III and IV consider geometry changes but 
retain only membrane forces. 
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I. Venkatraman and Hodge Solution [11] 



a. Clamped Edges 

The maximum deflection rate, at the plate center, is 
expressed in the form 



where G is an arbitrary function of time and 
n ,h.2+l/n 



)J= 2A 



(^) 



(B.I.l) 



(B.I.2) 



2n+l '2' 

The term X is obtained from the assumed form of Norton's law 
given as 

(B.I.3) 



® = <r> • 



When compared with the form given by equation A.l which was: 

(B.I.4) 



• Tr n 

e = Ko 



we have 



X = i 1/n 



(B.I.5) 



Substituting B.I.2, B.I.3, and B.I.5 into B.I.l: 



A* ^,a,2n+l ,, n , 2^'^^^” (2n+l) , n 
Wq = w* G(j^) Kaq ( ] 



(B.I.6) 



or 



w 



Ji^ 



n 



,a. 2n+l 



,V ,2^'^^/”(2n+l) ,n 

w* 241^^ ^ 



(B.I.7) 



which interprets B.I.l in terms of parameters introduced in 
Chapter III. 
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For comparison with the results in this paper we take 
G=l. The term w* is plotted as a function of n in reference 
[ 11 ] . 



b. Simple Support 

For simply supported edges the deflection rate at the 

plate center is of the form: 

9 2 

w = w* Ga (-35^ — (B.I.8) 
o 6u 

Substituting and rearranging as above this becomes: 

• 

w 

( 2 ) 

^£3 = W* G[-^ (B.I.9) 

-a. 2n+l 

where as before w* is plotted in reference [11] . 



82 



II. Odqvist Bending Solution [12] 
This result is of the form 



w = a n V (o) 
o 



(B.II. 1) 



where 



n+1 



„ _ 2n+l ,3a . 2n qa 

n '' . 2 ^ ha 

h c 



(B.II. 2) 



and V(o) is plotted as a function of n in reference [12] for 
the clamped and simply supported cases. The term is from 
the assumed constitutive relation: 



T a - s . . 

= l(_e.n-l 

eij 2'a ^ a 

c c 

which for uniaxial stress reduces to 



(B.II. 3) 



e = — 



n 



n 



(B.II. 4) 



and when compared with our notation: 
e = Ka” 



(B.II. 5) 



gives : 






(B.II. 6) 



Substituting these expressions into B.II.l and rearranging we 
have : 



w 



(• 



Kaq 



n 



-) 



(iilli)n ( 3 ) 2 V(o) 



-a. 2n+l 



n 



(B.II. 7) 
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Odqvist Membrane Solution [16] 



In this analysis the transverse deflection at the center 

of the plate is expressed in the form: 

2a h W(o) 

^ = — 2 (B.III.l) 

o q 



where 



W (o) 




P 

1 






2a h 
c 



U+-^c+[ 



(4n^+ 27n + 9) , ^21 
18(n+3) i ^ 



s 

8(n+3)s^'^^ 



s = [ ^ + 8(n+2)(n+3) ^] 



(B.III.2) 

(B.III.3) 

(B.III.4) 



(B.III.5) 



h is half the plate thickness, and C is tabulated as a function 
of n in the reference. In this form equation B.III.l gives the 
center deflection as a function of time. For the deflection 
rate we take the derivative with respect to time; 

2a h 5 

w = — ^ W(o). (B.III.6) 

o q 



Differentiating B.III.2 with respect to time; 



W(o) = - 



1 s 



u+ ic f , (B.iii.,) 



f rom B.III.4; 



s = -C 



8 (n+3) s 



,n+3 



(B.III.8) 
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Substituting B.III.5 and B.III .8 into B.III.7: 



00 / (4n^ +27n + 9)C^i 

W(o) = 2C(n+3) [l+tC+ 1 8 ( tit 3) 

[-^2 8(n+2)(n+3) — ] n+2 



(B.III.9) 



Substituting B.III.5 into B.III.2: 



P j 1 Ct ^ 

W(o) = -4 -[ - - ^To + 8(n+2)(n+3) 5=±_]n+2 



4 ' n+2 

s 

o 



rl+4c+ +27n +9) 2 , 

[I+ 7 C+ i 8 (n+ 3 ) ^ J 



(B. III. 10) 



Solving for t: 



t = 



r ( 4 W(o) >n+2 

8C(n+2) (n+3) ^ ' 2 B ' 

^ 1 



n+ 2 ' 



(B. III. 11) 



where B = 1 + 

18(n+3) 



(B. III. 12) 



Substituting 
W(o) = 



for t in B.III.9 
2C (n+3) 

4^+1 



P 



2 n +2 



(w(o) ) 



n +1 



p 2n+2 
1 

(W(o) 



where 

h 



2C B^~^^ (n+3) 
^n+1 



(B. III. 13) 



B. III. 14) 



Note that 4> here is that of reference tl 6 ] and bears no 
relation to the quantity <3> as used in Chapter II. 
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Odqvist's constitutive relation was; 



"ij = S 

c 



^ o , s . . 
3 ^ e.n-*l i] 



(B.III.IS) 



which reduces in uniaxial stress to: 



e = a" 

c 






giving : 



o = (i)!/" 

c 'k 



(B. III. 17) 



Substituting B.III.l, B.III.3, B.III.6, B. III. 14, and B. III. 17 
into B. III. 13 and rearranging: 



w 



(- 



Kaq 



n 



-) 



.a . 2n+l 



, 2h v n+1 
'w 

o 



(B. III. 18) 



Odqvist uses the symbol h to denote half of the plate thickness 
rather than the full thickness. If we call the full thickness 
h, we have a form analagous to that of the other results: 



w 

( — 8 ~) 



Kaq 



n 



.a. 2n+l 
'h' 



= 4 > 



f o. n+1 
'h ' 



(B. III. 19) 



where 



r 1 j_ 4n^ + 27n + 9 _2ii^t2 , ,-i\ 

i = 2C[H-ic t CM (iH-3) 

.n+1 



and C is tabulated as a function of n in reference fl6J . 
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I 



r 



IV. Membrane Solution for a Circular Plate, Onat and Yiiksel [17] 



Onat and Yiiksel assume that an initially flat circular 
plate deforms into a spherical surface upon application of a 
lateral hydrostatic pressure. Consideration of this geometry 
shows that the middle surface strain rate tensors have prin- 
cipal directions which are directed along the sphere's 
parallel, meridian, and the normal to the surface. The first 
two strain rates are positive and the third is given by 




da 



* ^ 01 
atan ^ 



(B.IV.l) 



For a spherical membrane: 

_ _ pR _ p 

*^1 ^2 ^ 2h 2h sina (1+cosa) 

o 

Creep deformations are governed by the relation: 



(B.IV.2) 



a tan | = f(a) (B.IV.3) 

where for Norton's law 

« 

f(a) = e (— )” = — a” = (B.IV.4) 

o a n 



The symbol p in B.IV.2 above represents the uniform pressure 
per unit area as used in [17] , is the diameter of the built- 
in support, and h^ is the original membrane thickness. 
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Figure B.IV.l 



From the geometry of the above figure we see 

w = R(1 - cos a) 
o 



D 

O T, 

= R sin a 



Equating R 



sin g 
(1-cosa) 



D 

o 



2 w 



o 



ot 

By trigonometric identity; tan ^ 

2w 

therefore, tan ^ =g— 



1-cosa 
sina ' 



Differentiating with respect to time 



d 

dt 



(tan^) 



a 

2 



2a 

sec 2 = 




or 



a 




D sec 
o 



2 « 

2 



(B.IV.5) 



(B.IV.6) 

(B.IV.7) 



(B.IV.8) 



(B.IV.9) 
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i 



I 



Substituting B.IV.2, B.IV.4, B.IV.8, and B.IV.9 into B.IV.3 



4w 



2w 



( 



:) (- 



r, 2d D 

D sec -- o 
o 2 



-) = K ( 



D p 
o 



2h 



(- 



n 



sina (1+cosa) 



and rearranging 
w 



(■ 



Kap 



n 



-) 



,a. n+1 
o 



w 

(— ) 
o 






1 .n 20 ti 

■ r) sec ■=•] 

sxna (1+cosa) 2 



(B.IV.IO) 



where 



D 

a = o 



The non-dimensional parameter resulting from the energy 
analysis is 
w 

( 



Kap 



n 



.a. 2n+l 
o 



Transforming B.IV.IO to this form: 



w 

Kap 



,a. 2n+l 



,a. n w 

^ ik 

o 



[i(— ^ r)^ sec^^] (B.iv.ll) 

sina (1+cosa) 2 



w 



O 3. 

Given and angle a can be found from B.IV.7 and values of 



w 



Kap 



n 



(-) 



2 x{-\-± computed as functions of ^ and 

o 



w 

_o 

1 

o 
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APPENDIX C 



Mathematical Details for Analysis with Combined Bending and 
Membrane Profile for the Clamped Circular Plate 

Experiments indicate that a clamped circular plate 
undergoing creep deformation from a transverse hydrostatic 
pressure assumes a velocity profile described by; 

w = w a-\-7) (c.i) 

O a. 

W W 

when is relatively small. At larger values of — the 
profile tends to the membrane form given by: 

w = w (1-fTr) (C.2) 

O a. 

These two forms can be combined by writing the velocity 
profile as: 

w = w^[aj^(l-|4)2 + a 2 (l-|t)I (C.3) 

where : 

tti + = 1 (C.4) 

and : 

0 < < 1 

0<a2fl (C.5) 

Equation C.3 reduces to the small deflection profile, C.I, for 
= 1 , = 0; and to the membrane profile, C.2, for = 0 , 

«2 = Using this profile and proceeding as in Chapter II: 
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(C.6) 



“r = “ “o'“l ^ ^ “2 i-1 



• w 

0 r 



• r 4 ,, r V . 2 , 

w [a, '^(1“^) + 

Ola a 2 a 



and 



K + 

r ^ e 



w 2 

^ (2a, (3-7^^) + 3a.] 
a. 1 a z 



6 ^ r 



w 2 

^ [2aj^(3-5|r) + 3a2l 



(C.7) 



(C.8) 

(C.9) 



In these equations and those to follow we may, of course, 
eliminate either or a 2 by using C.4. This was not done 
here, however, in order to keep the terms from the C.l profile 
and those from the C.2 profile directly visible and separated 
in each equation. 

For the radial deflection rate, u, we use the form: 
u = r (a-r) (Cj^+C 2 r) (C.IO) 



The extensional deflection rates are; 



e^ = u' + w'w' 



• u i 



Using as before: 

I w,w 

X = 



> C, = c* 



o o 

3 — 



w w ^ 
, o o 

' "^2 “ ^2 “ a ^ 



(C.ll) 



W W 

e = [ d-2x)c' + 

r a^ 1 



+ 16x^a^^a2 (1-x^ ) 



w w 

t(i-x)c' t 



(2x-3x^)c 2 + 16x^a^ (1-x^ ) ^ 

+ 4x^a2l 

(x-x2)c^l 



(C.12) 

(C.13) 
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These reduce to the appropriate expressions for the small 
deflection profile when = 1, = 0, and to the appropriate 



membrane profile expressions for 

w w 

+ 16 x^a^a 2 ( 1 -x^ ) 
+ 8x^a^ (1-x^ ) ^ + 



= 0 , and = 1. 

+ (■^x- 7X^C2 + 16x^a^ (1-x^ ) ^ 

+ 4x2a2] (C.14) 

+ (2x - |x^)c^ 

8x^aj^tt2 (1-x^ ) + 2x^a2](C.15) 



The general expression for moment per unit length derived 
earlier is; 



n 



„ ,4 ^ 1/n n 1 r^i ,bh|l+l/n 

” = <3K> H+T b 'sl*-"—! ■ T 2ntl)b 



a+ 



bh I 2+1/n 

2 Is 



. h I bh I 1+1/n . n 

* ll®-—! ^ l2H?i7b 



bh I 2+1/n, 



(C.16) 



and for extensional forces is: 



» = 1 (C.17) 



Where for M and N : 
r r 



a = 4j- + 7 e„ 

b = K + 

r 0 



and for and : 



• 1 * 

?e t T?r 

■=6 



Defining : 



Ac = 2aj^(3 - 7S^) + 3 u2 
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/ 



! 

i 



+ (|x 



Be = - ^x)c' ^ 



w 

7X )C2 



w 



+ 16x^a^ (1-x^ ) ^ ^ +16x^a^a2 (l~x^ ) ^ + 4x^0^ 

h 



w 

2^. 2 _0 
h 



D = 2a, (3 - 5x^) + 3a^ 
cl Z 



w 

~ 2 x)c| j:^ + ( 2 x - |x^) 



w w 

^'C2 jp + 8x^a^(l-x^)^ 



v; w 

+ 8x^a^a2 (l~x^ ) ^ + 2x^a^ ^ , 



v;e may write: 



M X = ^ hB-"!/" (^1 |B + iA 1 

r r 3K n+1 '■ 2A ' c ^ c‘ 

c 



7o".i:? ' \ ' I® + |2+Vn ^ , 5^1 B - 

(2n+l) ' c ^ c's ‘ 2A ' c ^ c 

c c 



1 + 1 /n 



n Sn^ - l) I Be - • t4aj(l-3x^)+2a2l(C.18) 



N e = “7T i [|b + ^A | 

r r 3K n+1 A ‘ c c‘ 



1 + 1 /n^ . rn- 9 x^r-' - 2 .+ ( 2 x - 3 x^)c,' — 



|B - iA 
c c 



J • [(l- 2 x)c 



1 h 



w 
_c 

'2 h 



w 



w 



+ 16x^a^ ( 1 -x^) ^ ~ + 16 aj^a 2 (l"X^)x^ 



w 

+ 4x^a| 



(C.19) 
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w 

.. /4 X 1/n n 1 . 2 + 1 /n / 1+1/n r 1 1 1+1/n 

= (^) h (^) 

c 



‘e"'9 ■ '3K^ 



n+1 



- T2H7I)^l ^ 2^1 i 



n 



P _ in |2+l/nn . 



(2n+l)D~^l ^ ^ ■*■ 202 ] (C. 20) 



»e-e = 4'^"" H?T 5^ 'l^c ^ 

c 



1 ^ , 1+1/n 



-|E, - 



w w 

] • [(l-x)c^ — + (x-x^)c^ ~].(C.21) 



The total energy dissipation rate is obtained by integrating 
these expressions over the deformed surface area of the plate. 
This gives: 



V = 2wa^ ,“o, 1+1/n 



'3K' 



n+1 






(C.22) 



where I is the integral of the bracketed expressions above 

w 

over X for a particular case (j~» n ) , and where the integral 
has been minimized with respect to and c^ for the value of 
and under consideration. 

Taking a small variation of w about its equlibrium 



w 



position for some particular value of 



dV » • 
dvv' o 



2tt ^ q 6w^dr 



(C.23) 



or 



• 1/n 2/4 . 1/n ,2+1/n ,1 > 1+1/n i s’ 



o 



2TTq6w^ ^ ( 1“^;^) ] rdr . 



2' a‘ 



(C.24) 
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Solving for w^; 



or 



and 



2+2/n 

' ' (i ) Vn ^ I 

3K 3-a^^ c 



,,a,2+l/n ,3Ka,l/n , ^"“l , 
"o ‘ "h> <Tn— > 

c 

in non-dimensional form: 

. '^o > _ 3 ,a.2n+l , ^ n 

V 4 42 I ’ 

Kaq C 









(C.25) 



(C.26) 



(C. 27) 
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APPENDIX D 



Mathematical Details for Analysis of a Simply Supported 

Circular Plate 



In calculating the elastic deflection of a simply 
supported circular plate subjected to a uniform lateral 
pressure Timoshenko and Woinowsky-Krieger [ 26 ] find that the 
deformation profile assumes the following form; 



w = w (1- —2-) (1 
o a 



1 + V 
5 + V 



(D.l) 



where w is the center deflection. If we make v = 1/2 and 
o 

take derivatives with respect to time, we get the deformation 
rate profile; 





W : 


= w (1 ~ ^ 
o 


) (1 - 


3 

11 




(D. 


.2) 


One 


may 


now compute 


: 












• 

K 

r 


w 

— “"W — — 7- 

a 


.28 

4 l 


36 

11 


r^l 


(D. 


. 3 ) 






. , w 

_ _w _ 0 

~ ~ a^ 


.28 

4 l 


12 

11 


1^1 ■ 


(D. 


. 4 ) 


and 


















r 


w 


[12 .. 
4 l 


42 

11 




(D. 


5 ) 




*^0 


W 

. 1 A _ 0 

+ 7 tC = — T 

^ r a ^ 


- 

‘11 


30 

11 




(D. 


6) 


Using, ( 


as in Chapte 


r II, 


the 


expression ; 








U : 


= r(a-r) (c^^ 


+ C2r) 






(D. 


7 ) 



for radial deflection rate where. 
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(D,8) 



c, = c* — 
1 la 



w w 
o o 



, o o 
= c! 



w w 

c 

TT 



'2 a 



and letting x = r/a: 



w w 

* * I ■ * I I o o 

e = u * + w* w = 



7 ^ [(l-2x)cj^ + (2 x-3x2)c^ 



+ x^ ( 



,28 12 



11 11 



x") 



. w w 

e,, = p = [(l-x)cj^ + (x-x^)c^] 



'0 



and 



Gj. + [(I - |x)cj^ + (|x - |x^)c^ 



(D.9) 

(D.IO) 



+x- 



(28 _ 12 2) 2i 

4l 11 ’ ^ 



(D.ll) 



w w 

Gq + [(I - 2x)cj^ + (2x - |x )c^ 







.28 12 
'll ■ 11 



x")2] 



(D. 12) 



These expressions may now be substituted into the equations for 
moment and extensional force resultants, (2.14) through (2.17). 
Defining ; 



A 



(il _ 42 2) 

'll 11 ^ ^ 



B 



= a - + (lx - - if=‘'>'ir 



E 



lU - ^ xM 



11 11 

(’^ - 2x)cj^ + (2x - fx^)c^ 



w 



w 2 Of, TO W 

_2 4 . ^ 12„2x 2 o 

2 h 2 'll 11 ’ h 
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( 



i 

‘i 




we may write: 



.. ♦ /4 ^ 1/n n 1 . 2 + 1 /n ,'^o^ 1+1/n r 1 i„ . i, 1 1+1/n 



n+1 



2A ' s ‘ s 
s 



n |B + + ^1 B - 

(2n+l)A ' s s's 2A ' s s ' 



, n |„ |2+l/n, r28 36 2 i 

+ 72STTT^I®s - > • ‘it - IT l 



(D.13) 



He = (If) 4i- (^l+l/n 1 (|b + |l+l/n 

r r 3K n+1 A ' s s ' 

s 



-|B, - • t (1 - 2x)cj^ + (2x - 3x^)c^ 



+ xM?? - ^x2)2 



11 11 



h 



(D.14) 



.« •, /4 ^ n 1 , 2 + 1 /n 1+1/n r 1 i t-. . ir, 

s 



1+1/n 



> 2 |e + io 1 2+1/n + i_ |e - 

(2n+l)D ' s s's 2D ' s s' 

s s 



+ |E_ - -iDjf I/"! • [ff - ^ xM 



n 



(2n+l)D| ' ^s ‘s's 



(D.15) 



• _ \ 1/n n , 2+1/n ,'^o. 1+1/n 1 



^0‘^O 



n+1 



h 






[ |e + ?d I 

D ' s s ' 






w 



w 



-|Es - • [(l--x)cj^ ^ + (x-x2)c^ j^] . 

(D.16) 
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< 



The internal energy dissipation rate is: 



V = 2Tra^ ^e'^e ^r®r 



(D.17) 



giving with (D.13) through (D.16): 



w 



V = 4r J 

3K n+1 a^ s 



(D.18) 



where is the integral of the bracketed geometric terms in 
the above expression. From the principle of virtual velocities 



— 6w = 2 ti /^q 6w rdr 

j ^ 0 

dw 

o 



(D.19) 



which gives: 



w 



1/n ,1 V 1+1/n _ 2/4 . 1/n ,2+1/n -r .. 

, (^) 2iTa^ (^) ' h ' 6w^ 



s o 



2-rrq 6w (1-fr-) (1“TT tdr = 2iTq6w ^ a^ . (D.20) 

On o. XX o. O 



So : 



= ( 



31 



2+2/n 



n 



A ,1/n . 2+1/n 22 
<3K> h 5 



(D.21) 



or w„ = I(|)2+Vn ,^,1/n g_l ,n 



o 



¥- 1 

5 s 



(D.22) 



or w 



Kaq 



o _ /3\ 2n+l 3 p 1 , n 



n 



4 "22 T 

5 s 



(D. 23) 
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APPENDIX E 



Illustration of Technique for Finding the Minimum Energy Integral 

As explained in Chapter II, the quantities and are 

determined by requiring that the internal energy dissipation 

rate (2.37) be minimized. Since the integration required to 

evaluate (2.37) must be done numerically, this minimization 

must be accomplished through a systematic search. The 

w 

minimization must be done for each case (p^f n) of interest; 
but, by considering a systematically varied series of cases, 
a set of curves such as those in Figures 2.2 and 2.4 can be 
obtained . 

The procedure used for minimization can be illustrated 

w 

by considering a typical case, n =» 5.0 and ~ 1*0 

clamped edges. The computer program listed in Appendix H, 

when provided an initial value for cj^ and a value for c^f 

will compute the integral, I^, for; the given c| and c^ values 

at regular intervals of 0.3. In the accompanying plots, the 

quantity $ rather than the integral I^ is plotted so that 

may be transferred directly to Figure 2.2. As = f (y~ ) t the 

■'‘c 

maximum value of corresponds to the minimum value of I^ . 

As seen in the plots, cj^ is varied until a peak is found for 

a particular value of c^- Peaks are then found for other 

values of c^ until a maximum is found with respect to c^. The 

resulting envelope of curves gives a maximum value of which 

gives a point on the plot of Figure 2.2. The procedure is 

w 

then repeated for the other values of ^ and n. 



100 




o 

iri 



o 



o 

fO 



q 

cJ 



q 



o 



lO o ^ 

to fO 



UJ 

o 



•e* 
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FIG. 



There is no mathematical assurance that the value of <J> 



found as above represents a unique and absolute maximum value. 
The initial area of search was that suggested for the analytical 
solution of the n = 1.0 case; and a very broad parameter search 
with respect to cj^ and c^ gave no evidence of any higher 
maxima. The agreement between the experimental results and 
the theoretical results suggests that the values of $ obtained 
are, in fact, the proper ones. 
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APPENDIX F 



Tabulation of Data from Plate Experiments 

The following tables list the data taken during the plate 
experiments discussed in Chapter I. A. All tests were taken at 
room temperature which generally was 75°F ±3°F with some points 
outside these limits. It is, of course, a well known fact that 
creep behavior in metals is affected by temperature. Within 
the relatively narrow temperature band here considered, however, 
such effects may be expected to be rather small and definitely 
secondary to the influence of pressure and geometry effects 
which were to be investigated in the plate tests. Examination 
of the data reveals no appreciable systematic influence of 
temperature in these tests. All times are in hours and minutes, 
all deflections are in inches, all pressures in pounds per 
square inch, and all temperatures in °F. 
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8 LB. PLATE 



TEST 1 



3.0 psi 



Time 



Hrs-Mins 


II 

o 




r=l" 


r=l^" 


r= 2 " 


r= 2 ^" 


00 

01 


. 0052 
53 


.0049 


.0033 


.0021 


.0005 


.0001 


03 


54 


49 


33i 


21 


05 


OI 7 


15 


57 


53 


36 


22 


06 


Ol'l 


32 


58i 


55 


37 


24 


06i 


Oil 


50 


60 


57 


38^ 


24 


07 


02 


1-00 


61 


57 


39 


24 


07 


02 


1-30 


63 


59 


40 


25 


07^ 


02 


2-00 


64 


61 


41^ 


26 


09 


02 


2-42 


66 


63 


43 


27 


09 


02 


3-00 


67 


64 


43 


2l\ 


09 


02 


4-45 


71 


67 


61 


29 


10 


02 


5-20 


12\ 


68 


47 


29| 


10 


02 I 


6-00 


73 


68^ 


48 


30 


10 


02 I 


6-30 


74 


70 


48^ 


31 


10 


02 I 


25-00 


92 


88 


62 


39 


14 


07 


47-00 


.0110 


.0105 


75 


47 


17 


08 


52-00 


14 


08 


78 


48 


18 


08| 


73-00 


23 


20 


86 


53 


20 


10 


100-30 


30 


26^ 


90 


55i 


21 


10 


125-40 


34 


30 


93 


57 


2li 


11 









TEST 


2 








8 LB. 


PLATE 




3. 


0 psi 




Time 


r=0 


II 

U 


r=l" 


r=l^" 


r=2" 


r=2|" 


00 


.0052 












01 


53 


.0045 


.0036 . 


0020 


.0007 


.0002 


03 


54 


45 


37 


21 


07 


02 


15 


56 


51 


38| 


23 


08| 


02 


32 


58 


53 


40 


23| 


09 


02 


50 


59 


53 


41 


24 


09 


02 


1-00 


59| 


53| 


41 


24 


09 


02 


3-50 


65 


62 


45 


28 


11 


02 


4-30 


67 


64 


46 


28 


111 


02^ 


5-00 


68 


64^ 


47 


29 


12 


03 


5-32 


69 


65 


48 


29 


12 


03 



Temp 

op 

72 



72 



72 

72 



72 



72 

72 

74 

76 

76 

76 

76 

76 



Temp 

op 



77 



76 

76 

76 

76 

76 

76 
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8 LB. PLATE 



TEST 3 



4.1 psi 



Time 




r=^" 




r=l^" 




r=2t" 


Temp 


Hrs-.Mins 


II 

o 


r=l" 


r=2" 


°F 


00 


.0352 














01 


.0388 


.0353 


.0280 


.0181 


.0060 


.0023 


69 


03 


.0405 


70 


93 


89 


63 


25 




07 


26 


85 


.0304 


97 


65 


26 




11 


33 


94 


11 


.0201 


67 


26 




18 


44 


. 0404 


19 


05 


69 


28 




30 


56 


15 


28 


10 


71 


29 


68^ 


47 


68 


26 


36 


16 


73 


29 




1-00 


75 


33 


41 


19 


74 


30 


00 


1-26 


84 


41 


47 


23 


75 


31 




1-49 


93 


48 


54 


27 


77 


31 


69 


2-30 


.0510 


64 


66 


34 


80 


36 




3-00 


17 


70 


72 


37 


81 


36 


69 


4-10 


32 


84 


81 


44 


83 


37 




5-00 


40 


91 


88 


48 


85 


37 


69 


5-30 


46 


98 


92 


51 


86 


. 38 




5-50 


49 


.0500 


94 


52 


87 


38 


69^ 



TEST 4 



8 LB. PLATE 



Time 


r=0 


r=r' , 


, r=l" 


00 


.0310 






01 


.0355 


.0311 


.0251 


03 


63 


25 


62 


06 


78 


38 


77 


11 


91 


50 


83 


15 


98 


56 


88 


22 


.0407 


64 


94 


30 


18 


73 


.0300 


45 


30 


85 


09 


1-00 


37 


91 


15 


1-41 


56 


.0410 


28 


2-00 


62 


15 


32 


2-30 


71 


24 


39 


3-00 


79 


31 


45 


3-33 


87 


39 


50 


4-06 


94 


45 


56 


4-30 


.0500 


50 


60 


5-00 


05 


55 


64 



4.1 psi 



r=l^" 


ii 


r=2V 


Temp 








op 








70 


0145 


.0055 


.0017 




51 


58 


18 




57 


60 


19 




61 


62 


21 




64 


63 


21 


70 


67 


64 


21 




72 


65 


22 




77 


68 


22 




80 


69 


23 


70 


88 


72 


24 




91 


74 


25 


70 


95 


76 


25 




97 


77 


25 


lOk 


0201 


78 


25 




04 


80 


26 


71 


07 


80 


27 




09 


82 


27 


71 
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8 LB 



PLATE 



TEST 5 



6.1 psi 



Time 

Hrs-Mins 


II 

o 


II 


11 


r=l^" 


r=2" 


II 


Temp 

op 


00 

01 


.0824 

.0869 


.0805 


.0614 


.0376 


.0150 


.0055 


79 


03 


96 


28 


32 


87 


55 


57 




06 


.0916 


47 


46 


95 


59 


59 




15 


44 


75 


66 


.0408 


64 


62 




25 


61 


94 


81 


17 


68 


64 




31 


68 


.0902 


86 


21 


70 


64 


79 


45 


87 


15 


98 


27 


74 


66 




1-00 


99 


28 


.0708 


34 


77 


66 


80 


1-30 


.1020 


48 


23 


44 


80 


68 




2-00 


35 


64 


34 


52 


84 


70 


80 


3-00 


59 


89 


53 


64 


90 


72 




3-36 


72 


.1000 


62 


69 


93 


73 


81 


4-45 


93 


20 


76 


79 


98 


75 




5-15 


.1107 


27 


86 


83 


.0200 


82 


81 


6-00 


19 


40 


93 


88 


03 


82 











TEST 


6 










8 LB. 


PLATE 






6.1 psi 






Time 


It 

o 




It 


r=l^" 


r=2" 


r=2i-" 


Temp. 












op 


00 


.0840 














01 


.0883 


.0842 


.0631 


.0379 


.0156 


.0073 


76^ 


03 


. 0908 


65 


49 


90 


61 


75 




07 


33 


86 


66 


.0400 


65 


78 




15 


55 


.0910 


83 


11 


70 


80 


76j- 


34 


85 


36 


. 0707 


24 


76 


85 


48 


98 


49 


17 


30 


79 


86 




1-00 


. 1006 


57 


22 


34 


80 


87 


76^ 


1-25 


21 


70 


35 


41 


83 


89 




1-50 


33 


81 


44 


47 


87 


92 


77 


2-00 


38 


87 


48 


49 


88 


93 




3-15 


61 


.1009 


68 


59 


92 


95 


78 


4-20 


78 


23 


81 


66 


96 


97 




5-50 


98 


43 


98 


76 


.0201 


.0100 


78 
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8 LB. PLATE 



TEST 7 



8.01 psi 



Time 

Hrs-Mins 


II 

o 


__ 1 •< 
r=-z 


r=l" 


r=li" 


r=2” 


r=2i" 


Temp 

op 


00 

01 


.1270 

.1313 


.1251 


.0943 


.0613 


.0244 


.0120 


73 


03 


39 


76 


65 


28 


50 


24 




06 


65 


98 


82 


39 


55 


29 




12 


88 


.1321 


. 1002 


52 


60 


32 


73 


34 


.1429 


57 


36 


73 


70 


38 




1-00 


54 


84 


57 


88 


76 


42 




1-27 


75 


.1403 


74 


98 


81 


46 


73 


2-00 


91 


22 


88 


.0708 


86 


48 




3-23 


.1524 


49 


.1114 


26 


93 


54 


73 


4-00 


34 


58 


22 


32 


96 


55 




5-15 


52 


76 


36 


42 


.0301 


58 


74 









TEST 


8 • 










8 LB. 


PLATE 




8 


.01 psi 






Time 


o 

II 

u 


r=i" 


1— 1 
II 

u 


r=li" 


r=2" 


r=2V 


Temp 
















op 


00 


.1252 












76 


01 


.1301 


.1247 


.0921 . 


0586 


.0237 


.0095 




03 


28 


74 


42 


99 


43 


98 




05 


43 


88 


55 . 


0607 


47 


99 




18 


89 


.1333 


91 


33 


57 


.0104 




30 


.1413 


54 


.1012 


45 


62 


10 


76 


45 


31 


70 


28 


55 


66 


12 




1-00 


45 


83 


39 


63 


69 


14 


76 


1-15 


56 


93 


48 


69 


72 


15 




1-30 


65 


.1403 


56 


75 


75 


16 




2-20 


89 


23 


76 


87 


80 


19 


76 


3-06 


.1504 


40 


90 


97 


85 


21 




3-35 


17 


51 


98 . 


0703 


88 


24 


75 


4-00 


24 


57 


.1104 


07 


89 


24 


76 


4-45 


36 


68 


14 


15 


93 


27 


76 


5-45 


50 


83 


26 


23 


96 


30 




6-00 


53 


86 


29 


25 


97 


30 


76 
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6 LB. PLATE 



TEST 9 



2.39 psi 



Time 

Hrs-Mins 


r=0 


r=4" 


r=l" 


r=l4" 


r=2" 


r=24’’ 


Temp 

op 


00 

01 


.0099 

.0100 


.0091 


.0067 


.0025 


.0012 


.0001 


74 


09 


124 


.01024 


744 


31 


134 


0l4 




15 


15^ 


05 


77 


32 


14 


Oil 




21 


18 


074 


784 


33 


144 


02 




30 


20 


11 


80 


344 


15 


02 




45 


23 


134 


824 


36 


16 


03 


74 


1-00 


254 


14 


84 


374 


16| 


03 




1-15 


274 


17 


86 


38 


164 


03 




1-30 


294 


194 


874 


384 


17 


034 




2-00 


324 


234 


90 


40| 


18 


05 


74 


3-17 


37 


30 


934 


43 


19 


05 




4-00 


384 


3 14 


95 


444 


19 


054 


74 


4-30 


394 


32 


964 


444 


19 


054 




4-45 


394 


32 


964 


444 


19 


054 




5-00 


4l4 


334 


97 


454 


194 


054 


74 


5-15 


4l4 


33| 


974 


4 54 


194 


054 




5-30 


4l4 


334 


974 


454 


194 


054 




6-15 


414 


364 


99 


47 


20 


054 


74 



Time 


6 LB. 
r=0 


PLATE 

r=4" 


TEST 

r=l" 


10 

2. 

r=l4 " 


39 psi 
r=2" 


r=24" 


Temp 

op 


00 

01 


.0097 
.0100 . 


0089 . 


0068 • 


0026 


.0010 


.0000 


74 


09 


09 . 


0100 


75 


31 


12 


00 


74 


15 


12 


03 


774 


32 


12 


oo4 




21 


14 


05 


79 


33 


13 


01 




30 


164 


07 


804 

824 


34 


13 


01 




46 


194 


11 


36 


134 


01 


74 


1-00 


214 


13 


84 


37 


14 


01 




1-15 


23 


134 


85 


38 


144 


01 




1-30 


244 


144 


86 


38 


144 


01 




2-00 


274 


18 


00 

00 


39| 


15 


02 


74 


3-15 


33 


23 


42 


16 


024 




4-00 


344 


254 


94 


434 


17 


024 




4-30 


35 


27 


944 


44 


17 


024 




5-00 


354 


28 


954 


44| 


17 


024 




5-15 


354 


28 


954 


444 


17 


024 




5-30 


354 


284 


96 


45 


17 


03 




6-15 


354 


294 


97 


454 


17 


03 
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6 LB. PLATE 



TEST 11 



4.1 psi 



Time 

Hrs-Mins r=0 r='^" 



00 


.0510 




01 


57 


.0498 


06 


79 


27 


12 


94 


42 


15 


.0600 


47 


21 


09 


54 


30 


19 


63 


45 


31 


74 


1-00 


40 


82 


1-15 




88 


1-30 


52 


93 


2-00 


62 


03 


2-30 


71 


11 


3-00 


79 


20 


4-02 


90| 


32 


4-30 


97| 


36 


4-45 


99 


39 


5-00 


.0701 


41 


5-15 


03i 


43 





6 LB. 


PLATE 


Time 


Hrs-Mins 


r=0 


r=i" 


00 


.0530 




01 


.0574 . 


0519 


06 


12 


52 


12 


26 


67 


15 


33 


71 


24 


43 


81 


30 


52 


89 


53 


72 . 


0607 


2-22 


.0708 


41 


2-45 


15 


47 


3-00 


21 


51 


3-31 


29 


59 


4-30 


40i 


70 


5-00 


44 




5-30 


51 


79| 


6-00 


54 


83 


6-12 


57 


85 



Temp 



r=li" 


r=li" 


CM 

II 

U 


op 








76 


0350 


.0265 


.0081 




71 


80 


86 




81 


87 


89 




85 


90 


90 


76 


91 


94 


92 




98 


99 


94 


76 


07 


05 


96 




13 


09 


91\ 


76 


18 


13 


00 




22 


15 


00^ 




29 


20 


02| 


76 


35 


24 


04 




41 


28 


06 


76 


49 


34 


08i 




53 


37 


09 


76 


55 


38 


10 




57 


40 


12 




59 


4li 


12 


76 



TEST 12 






r=l" 


r=li" 


4.1 psi 
r=2" 


Temp 

°F 


0399 


.0247 


.0092 


76 


24 


63 


98 




34 


71 


01 


76 


37 


73 


02 




44 


79 


05 




50 


82 


06 




65 


92 


10 


76 


94 


.0310 


17 


76 


99 


13 


18 




0502 


15 


19 


76 


09 


20 


21 




18 


26 


23 




23 


28 


24 


76 


27 


31 


25 




30 


34 


26 


76 


31 


34 


26i 
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6 LB . PLATE 



TEST 13 



6.1 psi 



Time 



Hrs-Mins 


r=0 


r=i" 


r=l” 


00 


.0880 






01 


.0927 


.0852 


.0650 


03 


51 


74 


68 


06 


68 


90 


80 


11 


85 


.0905 


94 


15 


94 


15 


.0701 


30 


.1017 


38 


20 


40 


29 


49 


28 


45 


33 


54 


33 


1-00 


45 


64 


42 


1-17 


57 


78 


53 


1-31 


66 


84 


59 


2-00 


83 


.0100 


71 


2-30 


94 


11 


80 


3-00 


.1104 


22 


89 


4-15 


26 


43 


.0807 


4-45 


33 


51 


14 


5-30 


45 


60 


23 



r=l^" 


(i 


r=2h“ 


Temp 

op 


.0441 


.0175 


.0071 


76 


53 


80 


73 




62 


84 


75 




71 


88 


78 




76 


90 


79 




90 


96 


81 


76 


96 


98 


83 




99 


.0200 


83 




.0505 


03 


85 


76 


11 


05 


87 




15 


06 


89 


77 


24 


11 


91 




31 


13 


92 


77 


37 


15 


93 


77^ 


50 


21 


95 


77i 


54 


23 


96 


60 


26 


.0100 





TEST 14 



6 LB. PLATE 



Time 


li 

o 


r=^" 


r=l" 


00 


. 0878 






02 


.0936 


.0865 


.0643 


06 


62 


88 


65 


15 


87 


.0915 


85 


31 


.1016 


42 


. 0706 


45 


32 


56 


18 


1-00 


44 


68 


28 


1-23 


60 


83 


40 


1-30 


66 


88 


45 


2-00 


79 


98 


55 


2-30 


91 


.1013 


66 


2-45 


97 


17 


70 


3-16 


.1106 


26 


77 


4-00 


18 


37 


87 


5-00 


30 


48 


98 


5-30 


36 


56 


.0802 



6.1 psi 



II 

u 


r=2" 


r=2T” 


Temp 

op 


.0438 


.0152 


. 0078 


80 


51 


57 


82 




65 


62 


86 




79 


68 


97 


80 


88 


71 


98 




94 


74 


.0100 


80 


.0503 


78 


03 




06 


79 


04 




13 


81 


05 


80 


21 


84 


07 




24 


84 


07 




30 


86 


08 


80 


36 


89 


09 




44 


92 


11 




48 


94 


12 


80 
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6 LB . PLATE 



TEST 15 



8.0 psi 



Time 




r=^" 




r=l^" 




<N 

II 

u 


Temp 


Hrs-Mins 


r=0 


r=l" 


r=2" 


°F 


00 


.1308 












75 


01 


.1354 


.1265 


.1009 


.0679 


.0255 


.0114 




03 


78 


88 


28 


93 


61 


16 




06 


98 


.1308 


44 


.0705 


66 


19 




15 


.1427 


36 


72 


25 


74 


23 


75 


24 


57 


54 


84 


36 


80 


27 




30 


63 


61 


92 


42 


83 


29 


75 


45 


73 


78 


.1108 


54 


88 


31 




1-00 


87 


91 


19 


63 


92 


33 


75 


1-15 


97 


.1403 


29 


70 


96 


35 




1-36 


.1510 


14 


40 


79 


.0300 


36 


75l 


2-00 


24 


27 


51 


88 


05 


38 




2-45 


44 


51 


69 


.0803 


11 


44 


76 


3-00 


50 


59 


74 


07 


12 


44 




4-20 


74 


81 


97 


25 


20 


47 


76 


5-04 


87 


95 


.1206 


33 


24 


48 




5-30 


93 


.1500 


12 


38 


26 


49 


767 


6-00 


.1600 


06 


18 


43 


28 


50 











TEST 


16 










6 LB. 


PLATE 






8.0 psi 






Time 


r=0 


II 

U 


r=l" 


r=li" 


r=2" 


r=2i" 


Temp 
















°F 


00 


.1295 












76 


01 


. 1342 


.1248 


.1009 


.0636 


.0265 


.0087 




03 


68 


78 


30 


53 


72 


89 




06 


88 


97 


47 


65 


77 


92 




15 


. 1417 


.1322 


73 


83 


86 


96 


76 


24 


34 


44 


88 


97 


91 


98 




45 


65 


68 


.1113 


.0716 


.0300 


.0103 




1-00 


78 


85 


25 


27 


05 


04 


76 


1-16 


90 


98 


37 


35 


09 


06 




1-30 


99 


.1406 


44 


42 


12 


07 




2-00 


.1517 


19 


60 


53 


18 


09 


76 


2-30 


31 


36 


71 


63 


23 


10 




3-04 


44 


47 


83 


73 


27 


12 




3-30 


53 


59 


92 


80 


30 


13 


76 
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TEST 17 



4 LB. PLATE 



0.82 psi 



Time 



Hrs-Mins 


II 

o 


II 


II 


00 


.0110 






01 


.0112 


.0106 


.0080 


03 


15 


08 


82 


06 


18 


13 


84 


15 


23 


16 


88 


30 


27 


21 


91 


1-00 


33 


26 


94 


1-34 


39 


32 


98 


3-00 


45 


38 


. 0103 


4-15 


53 


46 


08 


5-04 


56 


48 


10 


6-00 


57 


50 


11 



r=li" 




r=2^" 


Temp 


r=2" 


op 








75 


0050 


.0018 


.0003 




50 


18 


03 




52 


20 


04 




53 


20 


04 


74^ 


56 


21 


04 




58 


22 


04 


75 


59 


23 


05 




63 


25 


05 


74^ 


66 


26 


06 




67 


27 


06 


7 4? 


68 


27 


06 





TEST 



4 LB . PLATE 



Time 


r=0 


r=^" 


r=l" 


00 


.0127 






01 


.0131 


. 0122 


.0080 


03 


33 


25 


81 


06 


35 


27 


83 


15 


39 


31 


86 


30 


42 


33 


87 


1-00 


45 


37 


89 


1-30 


47 


39 


89 


2-00 


49 


42 


91 


2-18 


49 


42 


91 



18 


0.82 psi 






r=l^" 


II 


r=2V 


Temp 








op 








73 


0055 


.0020 


.0018 




56 


21 


19 




57 


23 


20 




59 


22 


20 


73 


60 


23 


21 




61 


24 


22 


73 


62 


24 


22 




63 


25 


22 




63 


25 


22 


73 
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TEST 19 





4 LB. 


PLATE 




Time 


Hrs-Mins 


II 

o 


11 

u 


r=l" 


00 

01 


. 0845 
.0871 . 


0820 . 


0653 


03 


83 


32 


64 


06 


92 


41 


72 


15 


.0906 


55 


83 


21 


11 


60 


88 


31 


18 


68 


93 


45 


24 


74 


98 


1-00 


29 


80 . 


0702 


1-31 


34 


86 


08 


2-00 


42 


92 


13 


2-30 


47 


98 


18 


3-00 


51 . 


0901 


22 



1.840 psi 

Temp 



r=li" 


II 

to 


r=2i-" 


op 

76^ 


. 0452 


.0170 


.0045 




55 


73 


46 




61 


75 


47 




68 


79 


48 




71 


81 


48 




76 


83 


49 


76^ 


80 


84 


49 




83 


86 


50 


76i 


88 


87 


50 




92 


89 


51 


76^ 


95 


91 


51 




98 


92 


52 


77 



TEST 20 





4 LB. PLATE 


Time 


r=0 r=|" r=l" 



00 


.0820 






01 


. 0848 


.0783 


. 0638 


03 


61 


98 


49 


06 


70 


.0804 


56 


12 


79 


14 


63 


18 


85 


24 


68 


26 


91 


31 


72 


30 


94 


33 


75 


1-00 


.0907 


46 


85 


1-15 


12 


49 


89 


1-30 


16 


56 


92 


2-00 


20 


58 


98 


2-30 


28 


67 


.0702 


3-30 


35 


76 


09 


4-19 


44 


84 


14 


4-30 


44 


86 


15 


5-30 


50 


92 


19 


6-00 


52 


94 


21 


6-33 


55 


96 


22 





1.840 psi 






r=li" 


r=2" 


<N 

11 


Temp 








op 








79 


.0438 


.0162 . 


0053 




46 


65 


55 




50 


68 


56 




56 


70 


57 


79 


61 


72 


58 




65 


74 


58 


79^ 


67 


75 


58 




74 


79 


61 


79i 


76 


80 


63 




79 


81 


63 


79 


84 


83 


64 




86 


85 


65 




92 


88 


66 


79 


96 


91 


68 




97 


92 


68 


79 


.0500 


94 


68 




02 


95 


68 


79 


03 


95 


68 
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TEST 21 





4 LB. 


PLATE 






4.1 psi 






Time 








r— 1 
11 
u 




r=2i" 


Temp 


Hrs-Mins 


II 

o 


II 


r=l" 


r=2" 


op 


00 


.1610 












79^ 


01 


.1651 


.1572 


.1333 


.0998 


.0377 


.0159 




03 


68 


87 


47 


.1010 


83 


61 




06 


81 


98 


57 


19 


88 


64 




12 


95 


14 


72 


31 


95 


67 




15 


.1700 


19 


77 


34 


97 


69 


79l 


25 


11 


33 


87 


44 


.0402 


70 




30 


14 


37 


90 


47 


04 


71 


80 


51 


27 


48 


.1403 


57 


10 


74 




1-00 


31 


54 


06 


61 


11 


75 


80 


1-30 


41 


63 


16 


68 


16 


79 




2-00 


49 


70 


23 


75 


19 


80 


80 


2-45 


58 


78 


30 


81 


24 


82 




3-00 


61 


81 


32 


83 


25 


83 


80 


3-30 


64 


86 


37 


88 


28 


87 


79^ 


4-03 


69 


88 


41 


91 


31 


89 




4-30 


72 


92 


45 


94 


32 


90 


79- 









TEST 


22 










4 LB. 


PLATE 






4.1 psi 






Time 


r=0 




r=l" 


r=li" 


r=2" 


r=2i" 


Temp 
















op 


00 


. 1580 












81 


01 


.1610 


.1538 


.1313 


.0973 


.0435 


.0136 




03 


30 


53 


29 


85 


42 


40 




06 


44 


68 


40 


95 


48 


42 




15 


63 


91 


56 


.1011 


56 


46 




24 


75 


.1602 


65 


19 


63 


48 




30 


80 


06 


69 


23 


65 


49 




48 


91 


18 


79 


32 


71 


51 




1-00 


96 


21 


83 


35 


73 


52 


81 


1-30 


.1706 


33 


92 


43 


78 


54 




2-00 


13 


43 


99 


50 


82 


55 


81 


3-05 


25 


54 


.1411 


59 


86 


59 




4-00 


33 


63 


17 


66 


90 


60 


81 


4-30 


37 


68 


20 


69 


95 


62 




5-05 


40 


70 


23 


71 


96 


62 


80^ 


6-00 


46 


76 


28 


75 


98 


63 


80 
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TEST 23 





4 LB. 


PLATE 




TIME 


Hrs-Mins 


II 

o 


(1 


ii 


00 

01 


.1875 

.1918 


.1820 


.1554 


03 


39 


35 


70 


06 


51 


47 


83 


15 


66 


66 


95 


25 


75 


78 


.1601 


37 


82 


85 


07 


45 


85 


89 


09 


1-00 


89 


94 


13 


1-15 


93 


94 


17 


1-33 


95 


.1901 


18 


2-00 


99 


02 


23 


3-00 


. 2009 


07 


27 


4-00 


14 


08 


30 


5-00 


17 


18 


39 


5-33 


20 


19 


39 



TEST 





4 LB. 


PLATE 




Time 


II 

o 


r=i" 


It 


00 


. 1925 






01 


.1970 


.1873 


. 1583 


03 


. 2000 


97 


. 1605 


06 


16 


.1908 


20 


08 


24 


14 


27 


15 


45 


37 


43 


26 


64 


60 


58 


30 


68 


62 


62 


45 


81 


72 


73 


1-00 


90 


81 


81 


1-30 


. 2103 


93 


93 


2-00 


14 


.2004 


.1702 


2-34 


24 


10 


10 


3-00 


29 


19 


15 


3-30 


35 


20 


20 


4-02 


40 


29 


26 


4-30 


45 


33 


29 


5-00 


49 


39 


34 


5-30 


54 


44 


37 



6.1 psi 



Temp 



r=l^" 


II 

to 


r=2i-" 


°F 

76 


1138 


.0488 


.0243 




53 


96 


46 




63 


.0503 


51 




77 


12 


54 




85 


18 


55 


76 


91 


21 


57 




93 


23 


58 




94 . 


26 


60 


76 


99 


27 


61 




1202 


29 


61 




05 


30 


62 


76 


09 


34 


63 


76^ 


12 


35 


64 


77 


17 


38 


71 




18 


39 


71 


77^ 



24 

r=li" 


6.1 psi 
r=2" 


r=2^" 


Temp 

op 


1187 


.0451 


.0212 


84 


1203 


62 


18 




14 


71 


22 




20 


76 


24 




33 


84 


28 


84 


47 


95 


32 




49 


96 


33 


84 


58 


.0502 


36 




73 


07 


38 


84 


73 


15 


41 




80 


20 


44 


84 


85 


25 


48 




89 


27 


49 


84 


94 


30 


50 




1300 


31 


51 


84 


03 


34 


52 




06 


36 


54 




09 


38 


55 


84 



115 



APPENDIX G 



Spectrochemical Analysis Report 



The reported analysis was made by Mr. Walter W. Correia 
of the STRNAD Spectrographic Laboratory at M.I.T. The sample 
labels, "4", ”6”, and ”8” are for the 4 lb. plate, 6 lb. plate, 
and 8 lb. plate respectively. The analysis samples were 
taken from tensile test specimens randomly selected from 
those actually used in the tensile experiments. From this 
data one may infer: 

4 
6 
8 



lb . 


plate 


Pb 


> 


99.67% 


lb. 


plate 


Pb 


> 


99.96% 


lb. 


plate 


Pb 


> 


99.29% 
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MASSACHUSETTS INSTITUTE OF TECHNOLOGY 
STRNAD SPECTROGRAPH 1C LABORATORY 
Rm. 13-4151 



Report of Spectrochemicol Analysis 

To. Peter T» Tarpgaard Charged To= 

Dote Received- 

Dote Reported-- 

Somple Description: £]? 



SAMPLE 


4 6 8 


SAMPLE 


4 6 8 


SAMPLE 


CO 

VO 


Ag 


VVFT 


VFT 


VFT 




Mg 


VFT' VFT 


VFT^ 




Ti 










Al 










Mn 










Tl 










As 










Mo 










U 










Au 










No 










V 










B 










Nb 










w 










Ba 










Ni 










Y 










Be 










Os 










Zn 










Bi 


T 


FT-T 


T 




p 


H 


H 


H 




Zr 










C 










Pb 










Ce 










Co 


WFT 


VFT 


WFT 




Pd 










Dy 










Cd 


T 




VFT 




Pt 










Er 










Co 










Rb 










Eu 










Cr 










Re 










Gd 










Cs 










Rh 










Ho 










Cu 


VVFT 


FT 


VFT 




Ru 










La 










Fe 










Sb 


FT-T 


FT-T 


T 




Lu 










Go 










Sc 










Nd 










Ge 










Se 










Pr 










Hf 










S i 


VFT 


VFT 


WFT 




Sm 










Hg 










Sn 


FT 




L 




Tb 










In 


T 








Sr 










Tm 










Ir 










To 










Yb 










K 










Te 




















Li 










Th 





















Remnrk*;: Sn In sample 6 approx« 



Vork by-. 
Reported 



by' 



Book No.?2«c8_|>oge No. -?.Q_ 



Key; VVFT <0.0001 % 

VFT 0.0001 % -0.001 % 
FT 0.00 1 %-O.OI% 

T 0.01 % -0.1 7o 
L 0.1 % - 1.0% 

M 1.0% -10.0% 

H >10.0 % 
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APPENDIX H 



Computer Programs 

I . Program for Clamped Edge Case 

The following listing is for the program used in the 
numerical evaluation of the energy dissipation rate integral 
for the clamped edge circular plate. The program consists of 
a main part and two subroutines. Subroutine SAMS evaluates 
the moment terms, ® radial position x, where 

X = r/a. Simlarly Subroutine JOES evaluates the extensional 
force terms, N^e^ + N^eg, at a position x. The main program 
uses five point Gauss quadrature to evaluate the total 
integral over x. The program is written in FORTRAN IV. 

The inputs to the program are n, c , c^/W^/h, and The 

input format is F10.5 which means each input quantity is 
alloted 10 spaces and if no decimal is inserted by the 
programmer the decimal is assumed to be to the left of the 
last five spaces in the field. The format for the input card 
is as follows: 

Field Quantity 



1-10 


n 


11-20 




21-30 


^2 


31-40 


w /h 
o' 


41-50 


a 



The last data card must have a negative number in the field 
1-10. This causes the program to exit. For each data card 
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the program computes and prints seven sets of output data 
v;here : 






(H.l) 



starting with the input value of c| plus .1. All other 
parameters remain constant through all seven iterations. 
The output is in ten columns as follows: 



Column 

1 

2 

3 

4 

5 

6 



8 

9 

10 



Quantity 

n 






c * 

w /h 
o' 



a. 



A , the integral of 
m ^ 

M K + 
r r 0 9 



A^, the integral of 

N e + 
r r 0 0 



A +A = I 
m n c 



. n 



where 4' = 



TTT 
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vO 
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z 
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> 






> 






> 






> 




to cn 






> 




> 






> 






> 
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9k 
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•k 


< 
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rH 
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rH 






rH 
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rH 






rH 






rH 


u 
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< 
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•k 
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2! 






y. 
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CD 








Z' 






z 






z 


o 


















UJ 






UJ 






UJ 






UJ 






UJ 




h- o 


r- 




UJ 




UJ 






UJ 






UJ 


o 
























9k 












» 






9k 




a * 


<r 




9k 




•k 






•k 






•k 


UJ 


< 
















r 






X 






X 






X 






X 




cr 


• 




X 




X 






X 






X 




f\j 


































'B: 










in no 


4- 






















Q 


u 




























9k 






9k 






9k 




oo 
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•k 






•k 






•k 


UJ 
























CM 






CM 






CM 






CM 




X ro 
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CM 






CM 


Q- 
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U 






U 






U 






U 




H- 00 


o 




U 




13 






U 






U 




rH 




o 




•> 




r-4 
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*. 






•k 






•k 






9k 




ft-i oo 
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16 EEA = ABS(EE) 

EEAl = EEA**( 1.0+1. 0/e^) 
EEA2 = (EEA^»*(2.0 + l.n/E.N) )* 
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II . Program for Simply Supported Case 



The following program for the simply supported case is 
exactly like that for the clamped case except that a different 
assumed velocity profile is used. 

The input data, in format F10.5, is as follows; 



Field 


Quantity 


1-10 


n 


11-20 




21-30 


"2 


31-40 


w^/h 



Here again the last data card must have a negative number in 
field 1-10 to cause the program to exit. The output data is 
exactly as in the clamped edge program discussed in Part A 
except that the column for is omitted and here 
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APPENDIX I 



Computation of Total Deflection 

The discussion in Chapter II was centered about determining 

the deformation rate, w^. In practical problems it is 

frequently the total deflection, w , at some time, t*, which is 

w 

of interest. It has been shown that w is a function of r^, 

o h 

and also that the material creep constants, K and n, may 

change with time. To compute w^, therefore, a piecewise 

linearization technique can be employed in which the defoannation 

process is broken into a sufficient number of linearized time 

segments to provide the accuracy desired. One such technique 

is discussed in this Appendix. 

In dividing the • deformation process into time segments 

some guidance can be obtained from the shape of the tensile 

‘ w 

test curVes and from the shape of the 4> versus ^ curves. 

Figures 2.2 and 2.4. The time segments should be relatively 

short during primary creep and may be longer when deformation 

is proceeding more slowly. As the deformation process 
w w 

proceeds, increases; and when ^ has changed enough to make 
a significant change in $ a new linear segment should be 
introduced. 

Having divided the total time of interest, from t=0 to 
t=t*, into n segments one can proceed as follows: 

Let t=t^ be the time at the beginning of the segment 
under consideration and t=t^ the time at the end. From 
tensile test data determine the appropriate values for 
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in the 



the creep constants, K and n, at the time t^, 
center of the given time segment, where; 



t 

c 




(I.l) 



Using these values of K and n; and using the value of 
w 

at the beginning of the time segment , compute 
with the appropriate equation from Chapter II, (2.53) for 
clamped edge or membrane or (2.57) for simple support. 
Multiply this value of w^ by the total time length of 
the segment, t^-t^, and add the result to the value of 
w^ at t^. This gives w^ at t^. The procedure is then 
repeated for the next segment, and the next, and so forth 
until t* is reached. 

In illustration, consider the case of Test 13, a 6 lb. 
plate under 6.1 psi pressure. Assume that it is desired to 
compute at t=8 hours. The total time is broken into two 
segments, a relatively short initial segment from t=0 to t=2 
hours, and a longer segment, where the deformation is more 
gradual, from t=2 hours to t=t*=8 hours. Values of K and n 
are determined at the median times of the segments , which in 
this case are t=l hours and t=5 hours. From the tensile test 
results in Chapter II the appropriate values are seen to be; 



t=l hour 


n=5 . 0 


K=4.1 X 


10-^^ 


(1.2) 


t=5 hours 


n=5 . 0 


K=1.7 X 


io-i« 


(1.3) 



In test 13 the initial deflection was .0880" giving: 
w 

^ = .840. (1.4) 

h 
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If test data is not available the initial deformation can be 
computed using elastic or plastic analysis. Using (1.2) and 
(1.4) and the appropriate value of f from Figure 2.2 in 
equation (2.48) the value of for the first segment is found 
to be .0090"/hr. Deflection at t=2 hours may now be determined 
(Wo)t=2 " -0880 + 2(.0090) = .1060 (1.5) 

which gives: 
w 



Using (1.3) and (1.6) with Figure 2.2 and equation 
find for the second segment that w^ = .0020"/hi^* 
at t=8 hours is: 



( 1 . 6 ) 

(2.48) we 
Deflection 



(w^)4._o = .1060 + 6(.0020) = .1180" (1.7) 

A plot of this two step linearization and of the measured 
experimental curve is shown in Figure I.l. Also shown is a 
three step linear approximation made with the technique 
described above with periods 0 £ t _< 1 hour, 1 £ t £ 3 hours, 
and 3 £ t £ 8 hours. Still closer approximations to the 
actual curve could be made by dividing the time period into 
more segments . 

In the practical application of these results it is 
important that initial deflection of the plate be known. 
Computing this initial deflection, however, may be no simple 
matter. Finite deflection analysis techniques for plastic and 
elasto-plastic deformation of plates are usually quite complex. 
For materials in which creep is important there is often a 
high degree of sensitivity to the speed of application of load 
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[28] which makes even the selection of appropriate elastic 
and plastic material constants a rather difficult task. In 
fact, such elasticity concepts as yield stress and Young's 
modulus may become rather meaningless for some materials under 
conditions where creep is important. 

In cases where the initial deflection is expected to be 
quite small and elastic, one could calculate initial deflection 
using the bending solution of Timoshenko and Woinowsky- 
Krieger [29]. For larger elastic deformation the large 
deflection solution of Timoshenko and Woinowsky-Krieger [30] 
is available. The large deflection elasto-plastic deformation 
of circular plates is treated by Ohashi and Murakami [31] [32] , 
and the large deflection of elastic-perf ectly plastic circular 
plates has been analyzed by Crose and Ang [33] . 
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